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Abstract 

We will provide detailed arguments showing that the set of Maxwell 
equations, and the corresponding wave equations, do not properly de- 
scribe the evolution of electromagnetic wave-fronts. We propose a non- 
linear corrected version that is proven to be far more appropriate for 
the modellization of electromagnetic phenomena. The suitability of 
this approach will soon be evident to the reader, through a sequence 
of astonishing congruences, making the model as elegant as Maxwell's, 
but with increased chances of development. Actually, the new set of 
equations will allow us to explain many open questions, and find links 
between electromagnetism and other theories that have been searched 
for a long time, or not even imagined. 



1 Short introduction 

The theory of electromagnetism, in the form conceived by J.C. Maxwell, can 
boast 130 years of honored service. It survived the severest tests, proving 
itself to be, for completeness and elegance, among the most solid theories. 
Very few would doubt its validity, to the extent that they may be more 
inclined to modify the point of view of other theories, rather than question 
the Maxwell equations. The trust in the model has been strong enough to 
obscure a certain number of "minor" incongruities and to incite the search 
for justifications in the development of other theories. 
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Nevertheless, even if the time-honored equations excellently solve com- 
plex problems, they are not able to simulate the simplest things. They are 
not capable for instance of describing what a solitary signal-packet is, one of 
the most elementary electromagnetic phenomena. Alternative models have 
been proposed with the aim of including solitons, but they did not succeed in 
gaining a long-lasting relevance, because they were based on deliberate ad- 
justments, that, accommodating specific aspects on one hand, were causing 
the model to lose general properties on the other. 

The development of modern field theory, which was very prosperous in 
the years 1930-1960, has magnified the role of the equations, giving them a 
universal validation in the relativistic framework. This progress came to a 
stop, leaving however the impression of being not too far from the goal of 
compenetrating electromagnetism and gravitation theory. 

We are going to make some statements that many readers will certainly 
consider heretic. We think that the various anomalies, which are present in 
the model, are not incidental, but consequences of a still insufficient theoret- 
ical description of electromagnetic phenomena. Actually, it is our opinion 
that the flaws are more severe than expected, and therefore, such a fun- 
damental "brick" of Physics needs extensive revision. The review process 
must be so deep that the entire setting necessitates re-planning from the 
beginning. On the other hand, if it were just a matter of small adaptations, 
this revision would have already been made a long time ago. 

We shall start to analyse some substantial facts, that at a practical level 
may be considered marginal, with the aim to evidentiate contradictions. We 
solve these problems by suitably redesigning the Maxwell equations. This 
will allow for the construction of a new model, solving all the inconsisten- 
cies and achieving the scope of a better understanding of electromagnetic 
phenomena. In a very natural way, the new approach also leaves the door 
more than open, to those links and generalizations that were expected to 
come from the Maxwell equations, but which, although vaguely insinuated, 
could never be realized in practice. 

None of the gracefulness that characterizes the Maxwell model will be 
lost. The reader who has the patience to follow our arguments through to 
the end, will discover that all the pieces find their exact place in a global 
scheme, with due elegance and harmony. We do not wish to say more in this 
short introduction. The model will be developed step by step, up to its final 
form, in order to let the reader appreciate the phases of its maturation. The 
mathematical tools used are classical, and maybe dated. On the other hand, 
our intention is to examine what would have happened to the evolution of 
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Physics, if our model was taken into consideration, in place of the Maxwell 
equations. We will elaborate and clarify many important concepts, leaving 
the path well clear for future developments, not considered here due to lack 
of time. 



2 Criticism of the theory of electromagnet ism 

In this section, we make some fine considerations regarding the evolution 
of electromagnetic waves, and the way they are modelled by the Maxwell 
equations. We start by pointing out deficiencies mainly at the level of math- 
ematical elegance. These will reveal other more severe incoherences. In the 
end, even taking into account the correctness, up to a certain degree of 
approximation, of the physical approach, our judgement will be rather neg- 
ative. As a matter of fact, in section 3, with the aim of finding a remedy to 
the problems that have emerged, substantial revision will be proposed. 

From now on, until section 11, we assume that we are in void three- 
dimensional space. As usual, the constant c indicates the speed of light. In 
this case, the classical Maxwell equations are: 

— = c 2 curlB (2.1) 

divE = (2.2) 

^ = -curlE (2.3) 

dt y ' 

divB = (2.4) 

where the vector field E is dimensionally equivalent to an acceleration mul- 
tiplied by a mass and divided by an electrical charge; while B is a frequency 
multiplied by a mass and divided by a charge. 

The above equations are supposed to be satisfied point-wise at any in- 
stant of time. Their solutions are assumed to be smooth enough to allow 
differential calculus. Therefore, discontinuous or singular solutions are not 
allowed. The equations (|2.2|) and (|2.4|) could be considered unnecessary, 
since they are easily deduced from (|2.1|) and (|2.3j) respectively, after apply- 
ing the divergence operator. Later on, for the reasons that we are going to 
explain, we will question the validity of 1)2.2(1 and (|2.4|) . As a consequence, 
the entire formulation will lose its credibility. 
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As far as the evolution of an electromagnetic plane wave (with infinite 
extent and linearly polarized) is concerned, we have no objections to make. 
In Cartesian coordinates, a monocromatic wave of this type, moving along 
the direction of the z-axis, is written as: 

E = {csmuj(t- z/c), 0,0) B = (0, eanu(t - z/c), 0) (2.5) 

In this case, the Maxwell equations are all satisfied point-wise. 

The next step is to examine the case of a spherical wave, which is far 
more delicate. The wave could be generated by an oscillating dipole of 
negligible size. However, the way the wave is produced and supplied is not 
of interest to us at the moment, being more concerned with analyzing the 
geometrical aspects of its evolution at a distance from the source. 

Let us denote by P = E x B the Poynting vector. It is customary to 
assume that E and B are orthogonal, and that the wave-front propagates 
at constant speed c, through spherical concentric surfaces. One may argue 
that perfect spherical waves do not exist in nature. Nevertheless, for the 
sake of simplicity, we maintain this hypothesis which can be removed later, 
without modifying the essence of our reasoning. 

We are basically confronted with two possibilities. In the first one, the 
Poynting vector follows exactly the radial direction. This means that E 
and B locally belong to the tangent plane to the wave-front. In such a 
circumstance, as detailed below, we are able to show that (|2.1|) and ()2.3|) 
cannot be both satisfied everywhere. More precisely, it is known that (|2.1|) 
and H2.3|) are true up to an error that decays quadratically with the distance 
from the source. Since the intensity of a spherical electromagnetic wave 
only decays linearly in amplitude, the above mentioned inaccuracy has no 
influence on practical applications. However, we record a first negative mark. 

The second possibility is that, in order to satisfy all the set of Maxwell 
equations, we loose the orthogonality of the Poynting vector with respect to 
the wave-front surface. This is a more unpleasant situation, considering that 
the Poynting vector represents the direction of propagation of the energy 
flow. The lack of orthogonality between the wave-front tangent plane and 
the direction of propagation violates the Huygens principle (recall that we 
are in vacuum) , leading to a deformation of the front itself. As we will check 
later, this results in relevant defects in the development of the wave-shape. 

Let us study the problem more in detail, by taking into account the 
transformation in spherical coordinates: 

(x, y, z) = (r sin^cos^, rsinc/>sin#, rcos0) (2-6) 
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with < 9 < 2ir, < do < it and r large enough. We look for vector fields 
having the following form: 

B = (0, 0, u) E = 0, w, 0) (2.7) 

where u, v, w are functions of the variables t, r and <f> (no dependency on 
9 is assumed). In (|2.7j) . the first component of the vectors is referred to 
the variable r, the second one to 4>, and the third one to 9. The unknowns 
in the system of Maxwell equations reduce from six to three. Choosing a 
more general form for the fields only complicates the computations, without 
adding anything to the substance. 

We start by observing that equation (|2.4|) is immediately satisfied. More- 
over: 

' u cos 6 us ,u 

— ^+ ' ~ - 
r sin <p r r 



/ u cos <p us ,u . \ fn 
curlB= — —7 + ^, -(- + u r ), (2. 
V r sin m r r > 



divE = v r H 1 f H ^ (2.9) 

r r sin <p r 

curlE = fo, 0, w r + --^) (2.10) 

Therefore, the equations in spherical coordinates become: 

U t = ~(wr + ™) + V ~f (2.H) 

c 2 / cos 



c / cos cp \ 

-(u ^4 + us) 2.12 

r V sinffl / 



v t 

r v sin 



'U 

w t = - C ' 

'> r 



(- + n r ) (2.13) 



To avoid discontinuities, we must introduce the following boundary con- 
straints: 

u(t,r,0) = u(t,r,vr) = w{t,r,0) = w(t,r,7r) = (2.14) 
dv dv 

^ (t ' r ' 0) = ^ ( *' r ' 7r) = ° (2 ' 15) 

In the case of the pure radiation field of an oscillating dipole, when r is 
sufficiently large, one usually sets v = and w = cm. Within this hypothesis, 
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the two equations Q2.11JI and (|2,13|) are equivalent. They bring us to the 
general solution: 

w{t,r,<j>) = cu(jt,r,</>) = ^f{4>) g(t-r/c) (2.16) 

where / (with /(0) = /(vr) = 0) and g are arbitrary functions (the only 
restrictions apply to their regularity). Among these solutions there is the 
one corresponding to f((j>) = sin</>, which is often present in classical texts 
(see for instance p. 284), being the one with more physical relevance. 
Nevertheless, we unfortunately note that equation (|2.12|) is compatible with 
v = only when: 

B = (0, 0, — g(t - r/c)) (2.17) 
V r sm <fi > 

E = (0, — °— g(t-r/c), o) (2.18) 
\ r sm (p ' 

which manifest singularities at the points corresponding to 4> = and <f> = ir. 
In general, we have the following statement: 

divE = \g{t- r/c) {^-m + f'(<t>)) =0 & f^) = -^J ( 2 - 19 ) 

Note that such a strong singularity at the poles cannot be removed only by 
requiring the wave-front not to be perfectly spherical. 

We observe that / can be taken in such a way that divE is vanishing at 
the poles (for example f(<fi) = (sin</>) 2 ), but not in the proximity of them. 
In addition, we observe that, if / is regular with /(0) = f(ir) = 0, for any 
fixed r, the points in which the divergence of E does not vanish belong to 
a bidimensional set whose measure is different from zero. For instance, if 
f{4>) = sin0, we find out that divE is proportional to cos(/>, so that this 
set consists of all points of the sphere of radius r, with the exception of 
the equator. It is certainly true that even if the divergence is not zero, it 
is negligible when designing, for instance, a device like an antenna. This 
argument, however, is not going to be valid here, since we would like to 
carry out an in depth analysis of what is really happening in the evolution 
of an electromagnetic wave, compared to what the Maxwell theory is able 
to predict. 

Let us now follow a different path and try to find other solutions, of the 
form given in (|2.7|) . satisfying the set of all Maxwell equations (including 
divE = 0). If we do not want / to be singular somewhere, we have to accept 
that v is different from zero. This means that E has a radial component, 
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so that the Poynting vector cannot be perfectly radial. We have to better 
check what happens in this last case. 

It is well-known that the the Maxwell equations lead to: 

1 (>>E AE and 1^5 = AB (2.20) 



c 2 dt 2 c 2 dt 2 

The above are usually called "wave equations", but, shortly, we will see that 
this name is not appropriate. The terminology is correct only if the fields 
involved are scalar. By deriving (|2.11|) with respect to time and using (|2,12|) 
and (|2.13|) . we arrive at the equation: 



\uu = ^(r\) r + ^ (^T ( " M,lo) " ) S2 ' 21) 



corresponding to the third component of the second equation in (|2.20)) in 
spherical coordinates. 

It is worthwhile noting that (|2.21|) is not the wave equation for the scalar 
field u in spherical coordinates, due to the fact that in this framework the 
Laplacian of a vector field is not the Laplacian of its coordinates (even if 
only one of them is different from zero). The wave equation for u reads as 
follows: 

-o U tt = -^(r 2 U r ) r + {US SU1 = Au (2.22) 

cr r A r A sin q> 

This is not a trivial warning, since many texts in electromagnetism erro- 
neously confuse (|2.22|) with (|2.21|) . Implicitly, we made the same mistake 
before, when looking for d'Alembert type solutions of the form ()2.16|) . gen- 
erating, for this reason, solutions not compatible with all the Maxwell equa- 
tions. 

By separation of variables, for any k > 1 and any n > 1, we discover 
that ()2.21|) admits the following basis of solutions: 



r 2 cos(cfct) J n+ i{kr) sin0 P^coscj)) 

r~2 sin(cfct) J n+ i{kr) sm<p P^(cos <j>) 
r^a cos{ckt) Y n+ i(kr) sin <j) P' n {cos 4>) 

s'm{ckt) Y n+ i{kr) sin0P^(cos0) (2.23) 

where J n+ i and Y n+ i are Bessel functions of first and second kind respec- 
tively, while P n is the n-th Legendre polynomial. 



7 



A classical reference for Bessel functions is jllj . It is important to note 
that the solutions given in at page 127, for the scalar wave equation in 
spherical coordinates, differ from the ones shown in ()2.23|) . The reason is 
that the functions in jTJ (having P n (cos0) in place of sincp P^(cos^)) are 
those solving (|2.22|) . which is not the vector version of the wave equation, 
as we already mentioned. 
For example, if n = 1 we have (see p. 54): 

t fu \ fsmkr \ 

J a ykr) = \ — — — : cos kr 

2 V irkr \ kr J 

Y i {kr) = {^k~r +sinkr ) Pi ' (cos0) = 1 (2 - 24) 

In order to understand what the solutions in (|2.24j) look like, it is standard 
to introduce some approximation. Thus, for n = 1 and r large, by taking 
the combination r -1 / 2 (sin(cfci) Js/z^kr) + cos (ckt)Y 3 / 2 (kr)) sin0, up to mul- 
tiplicative constants, it is possible to get asymptotically the monocromatic 
solution u = r^ 1 sin0 sin k(ct — r) (compare to (|2.16[) ). up to an error which 
decays quadratically with r. Once again, one ends up with something simi- 
lar to a travelling wave, although some cheating has been necessary (that is 
equivalent, in the end, to replacing once again (|2.22|) by (12.211) ). 

On the other hand, suppose that u is evaluated exactly as linear com- 
bination of the functions in (|2.23j) . Then, one recovers v e uu by (|2.12|) and 
()2.13|) . through time integration. Successively, it is possible to compute the 
Poynting vector: 

P - c 2 (u J(^ + u^jdt, ^ Jl^-^^ +u<fijdt, Oj (2.25) 

which has, as expected, a non radial component. Now, let us fix r and study 
the behavior, by varying (j), of the two components of P. In particular, we 
are interested to see what happens near the poles (4> = or eft = ir). We start 
by noting that, for any n > 1, the term P^(cos 4>) tends towards a finite limit 
for <f> -» or cf> -* it (recall that P^(±l) = \(±l) n+1 n(n + 1)). Therefore, 
according to (|2.23|) . the first component in (|2.25j) behaves as (sin(/>) 2 near 
the poles. It is a matter of using known properties of Legendre polynomials, 
in particular the differential equation: 

(sin0) 2 P^(cos0) - 2 cos <f> P' n (cos (j)) + n(n + l)P n (cos (j)) = (2.26) 

to check that the second component in (|2.25|) behaves as sin cp near the poles. 
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We are ready to draw some preliminary conclusions. Let us note that 
finally divE = 0, hence all the Maxwell equations are satisfied. As already 
remarked, it has been necessary to keep the nonradial component of P. 
Surprisingly, for any fixed r, such a nonradial component prevails on the 
radial one, when approaching the poles. This implies that the shape of 
the wave-fronts does not resemble a sphere, but rather a kind of doughnut 
with the central hole reduced to a single point. The parts of the fronts 
corresponding to the internal side of the doughnut, progressively stratify 
along the z-axis. We do not see a chance of recognising any sort of Hyugens 
principle here. This is not what we would call a travelling wave. It may 
be argued that this behavior is due to the influence of the source located 
at r = 0. But, if we stop the source, the wave-fronts already produced 
continue to develop. If their motion is ruled by the Huygens principle, the 
hole should fill up quickly, and each front should transform to something 
rounded which is almost a perfect sphere. The problem is that, during this 
smoothing process, the vector fields E and B are not compatible with both 
the constraints divE = and divB = 0. The clue is that "wave equations" 
in vector form have nothing to do with real waves. 

Some mild analogy between the Maxwell equations and the eikonal equa- 
tion, governing the movement of the fronts, was devised a long time ago. The 
equivalence is valid within the limits of geometrical optics (see [Hj, p. 110). 
In spite of this, examining the behavior of the fronts, our impression is that 
their natural evolution is in conflict with all restrictions imposed by Maxwell 
equations. This statement will be clearer as we proceed with our study. The 
right connections with the eikonal equation will be defined in section 10. 

The different situations analyzed up to now are summarized in figure 1, 
which should clarify our point of view: either we keep the singularity at the 
poles (manifested by infinite amplitude of the fields or strong geometrical 
distorsion), or we allow the divergence of the electric field to be different 
from zero. To sustain this proposition, let us collect other elements. 

Some confusion usually arises when one tries to simulate the evolution of 
a "fragment" of wave. We examine the case of the plane wave given in 1)2.5(1 . 
For any fixed z, we can cut out a region O in the plane determined by the 
variables x and y, and follow its evolution in time. For simplicity, O can be 
the square [0, 1] x [0, 1]. Inside SI we assume that the electromagnetic fields 
evolve following 1)2.5(1 . in full agreement with Maxwell equations. Outside f2, 
the fields E and B are supposed to vanish. The question is understanding 
what happens at the boundary dO, of O. It is not difficult to realize that, 
on the sides {0}x]0, 1[ and {l}x]0, 1[, curlB and divE become singular, 
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Figure 1: Qualitative behavior of the held E as a function of the angle <f>. Case 1: 
E = (0, r~ 1 (smcp)g(t — r/c),0), the wavefronts are perfect spheres, but divE ^ 0. 
Case 2: E = (0, (r sin (f>)~ 1 g(t — r/c),0) the condition divE = is satished, but 
there are singularities at <j> = and </> = 7r; Case 3: tie corresponding Poynting 
vector is given in 12.25$ . the divergence of the electric field is vanishing, but the 
wave-fronts are far from being spherical surfaces. 

producing concentrated distributions. Similarly, on the two sides ]0, l[x{0} 
and ]0, l[x{l}, the quantities curlE and divB present singularities. 

Some readers may complain because discontinuities of the fields may 
not exist in nature. Commonly, the right way to proceed is to consider a 
thin layer around dQ, where the solution given by ()2.5|) smoothly decays 
to zero. Then, one lets the width of the layer tend towards zero. This in 
general allows us to determine special relations to be satisfied on d£l (in place 
of the Maxwell equations, which are meaningless there). Unfortunately, the 
procedure presents some drawbacks. Let us first assume that the wave- fronts 
shift along the z-axis maintaining their squared shape. We also assume 
that the fields E and B are orthogonal and smoothly decaying to zero in a 
neighbourhood of <9$7 (like for instance in figure 2). Our conjecture is that 
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there exists at least one point where Maxwell equations are not all satisfied, 
because divE and divB cannot both be zero at the same time. Actually, 
examining figure 2, we discover that there are infinite points where either 
divE ^ or divB 7^ 0. These points form a set whose area is different from 
zero. We are free to try other configurations by modifying the orientation 
of the vector fields at each point near <9f2, but we always arrive at the same 
conclusion: some rule of Physics breaks down when approaching <9f2. Now, 
the question is: if we do not know what the governing rules are in the layer 
around <90, how can we go to the limit for the size of the layer tending to 
zero? 

Another possibility is that the wave- fronts, due to the strong variation 
of the fields near the boundary of O, are forced to bend a little. The electro- 
magnetic fields are no longer on a plane, so we could probably find out the 
way to enforce all the Maxwell equations. However, this implies that the 
Poynting vectors are not parallel to the z-axis anymore. Thus, the shape of 
is going to be further modified during the evolution. A little diffusion is 
bearable, yet our impression is that the wave-fronts would rapidly change 
their form. The more they bend, the faster they produce other distorsion. 
This is in contrast for instance with the fact that neat electromagnetic sig- 
nals, of arbitrary transversal shape, reach our instruments after travelling 
for years between galaxies. The only acceptable rule is that all the Poynting 
vectors must stay orthogonal to the fronts and parallel to the actual direc- 
tion of movement; if this does not happen the wave quickly deteriorates, 
fading completely. 

To prove what we claimed before, we show using very standard argu- 
ments that it is not possible to construct solutions to Maxwell equations, 
having finite energy and travelling unperturbed at constant speed along a 
straight-line. We assume that the speed is c and the straight-line is the 
z-axis. Without loss of generality, such a signal-packet is supposed to be of 
the following type: 

E = (E 1 (x,y),E 2 (x,y),E 3 (x,y)^ g(t-z/c) 

B= ( K B 1 (x,y),B 2 (x,y),B 3 (x,y)) g(t-z/c) (2.27) 

where g is a bounded function and all the components E\, E 2 , E3, B\, B 2 , 
B3 are zero outside a bidimensional set Q. It is not difficult to check that 
(|2.1jl and 1)2.3)1 only hold when E and B are identically zero. Actually, it 
is straightforward to discover that £3 e B3 must be constant (and the sole 
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Figure 2: Example of electromagnetic Geld smoothly reducing to zero at the bound- 
ary of the square. The Poynting vector is orthogonal to the page at each point. 
The Maxwell equations are satisfied in the central part. Instead, approaching the 
boundary, the divergence of the fields turns out to be different from zero. 

constant allowed is zero). Then, one finds out that E±, E2, B\, B2 must be 
harmonic functions in £1. Since they have to vanish at the boundary, they 
must vanish everywhere. 

Due to the above mentioned reasons, solitonic solutions are not described 
by the classical theory of electromagnetism. Efforts have been made in 
the past to generalize the Maxwell model, in a nonlinear way, in order to 
include solitons. Just to mention an example, the Born-Infeld theory (see 
[3]) predicts the existence of finite-energy soliton-like solutions (that have 
been successively called Blons). These last equations have no relation with 
the ones we are going to develop in this paper. However, they point out the 
necessity of looking for nonlinear versions of the model. We will come back 
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to the subject of solitary waves in section 5. 

In many applications, a standard approach is to reconstruct the bidi- 
mensional profile of the fields inside O with the help of a truncated Fourier 
series. This is accomplished by a complete orthogonal set of plane waves, 
each one carrying a suitable eigenfunction in the variables x and y. We 
must pay attention, however, to the fact that these eigenfuntions are of the 
periodic type. Therefore, they reproduce the same profile, not only inside Q, 
but in a lattice of infinite contiguous domains. In this way, the represented 
solution turns out to have infinite energy. Considering only one of these 
profiles, thus forcing to zero the solution outside ft, unavoidably brings us 
again to a violation of the Maxwell equations near the boundary of f2. Some 
clarifying comments on this issue can be found in [H], p. 42. 

We recognize that the techniques based on Fourier expansions provide 
excellent results in many practical circumstances, as for example the study 
of diffraction. Nevertheless, in this last case and in the ones treated before, 
it is necessary to adapt the solutions, introducing some approximation, if we 
want them to correspond to the real phenomenon. Indeed, these adjustments 
are within the so-called limits of the model. Hence, we could just stop our 
analysis here, with the trivial (well-known) conclusion that the Maxwell 
model is not perfect. We believe instead that the discrepancies pointed out 
are not just imperfections, but symptoms of a more profound pathology 
affecting the theory of classical electromagnetism. 

What we learn in these pages is that there are plenty of simple and 
interesting phenomenon, which are inadequately explained by the Maxwell 
model, because the equations impose too many restrictions. Consequently, 
the idea we shall follow in the next section is of weakening the equations, 
with the aim of widening the range of solutions. 



3 Modified Maxwell equations 

The demolition process is finished, now it is time to rebuild. To begin, we 
propose the following model: 



<9E 
~dt 



c 2 curlB 





dB 

~dt 



cur IE 



(3.2) 
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divB = 



(3.3) 



that will be further adjusted in the subsequent sections. The norm | • | is the 
usual one in R 3 , i.e.: \{x,y,z)\ = \J x 1 + y 2 + z 2 . We define J = P/|P| = 
(E xB)/|Ex B|. Note that, when P = 0, the direction of J is not determined 
(see also the comments at the end of section 7) . The vector J is supposed to 
be adimentional (or, equivalently, P/|P| is multiplied by a constant, equal 
to 1, whose dimension is the inverse of the dimension of P). Consequently, 
cJ is a velocity vector. 

As the reader may notice, the "awkward" relation divE = has been 
eliminated. It is also evident that in all the points in which divE = 0, we 
find again the classical Maxwell system. This states that the solutions of 
C3J-!E3-C3|l-lE3J are also solutions of (|3TT|) - (|33)) - (PP|) . Therefore, the 
replacement of (|2.1j) - (|2.2j) by (|3.1|) brings us to the property we wanted, that 
is the enlargement of the range of solutions. 

Afterwards, we have to understand and justify what is happening from 
the point of view of Physics. Let us recall that we are in empty space, 
and that there are no electrical charges or masses anywhere. We are only 
examining the behavior of waves. In spite of that, we pretend that there 
may be regions where divE ^ 0. The situation is not alarming, since we 
checked that the condition divE ^ is quite frequent in the study of waves. 
Anyway, such a hypothesis is acceptable, as long as it is coherent with the 
basic laws of Physics. First of all, we observe that the equation ()3.1|) has 
been obtained by adding a nonlinear term to (|2,1|) . The term has strong 
analogy with the corresponding one of classical electromagnetism, appearing 
on the right-hand side of 1)2- lj) as a consequence of the Ampere law, and due 
to the presence of moving charges. In fact, by setting p = divE, the vector 
pcJ can be assimilated, up to dimensional constants, to an electric current 
density. Thus, even if in our case there are no real charged particles, we have 
to deal with a continuous time-varying medium, consisting of infinitesimal 
electrical charges, living with the electromagnetic wave during its evolution. 
Moreover, the added term does not compromise the theoretical study of a 
functioning device like an antenna, since, at a certain distance, the quantity 
divE is negligible. 

By taking the divergence of (|3.1|) . we get a very important relation: 



which is, actually, the continuity equation for the density p = divE. The 
equation (|3.4|) testifies to the presence of a transport, at the speed of light, 



dp 



c div(pJ) 



(3.4) 
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along the direction determined by J. Hence, something is flowing together 
with the electromagnetic fields; something that later, in sections 9 and 10, 
will be compared to a true mechanical fluid. On the other hand, this was 
also the interpretation at the end of the 19th century, before the theory of 
fields was rigorously developed. The fluid changes in density, but preserves 
its quantity, as stated by the continuity equation. It is extremely signif- 
icant to remark that this property comes directly from (|3.1|) . so it is not 
an additional hypothesis. In section 12, based on the density p, we will 
construct a mass tensor that, due to QM.lj) . can be perfectly combined with 
the standard electromagnetic energy tensor. The skilful reader has already 
understood that this will allow us to find the link between electromagnetic 
and gravitational fields. 

We are now going to collect other properties about the new set of equa- 
tions. Considering that E x B is orthogonal to both E and B, a classical 
result is obtainable: 

1 8 

- — (|E| 2 + c 2 |B| 2 ) = c 2 (curlB-E-curlE-B) = - c 2 divP (3.5) 

where the quantity |E| 2 +c 2 |B| 2 , up to a multiplicative dimensional constant, 
is related to the energy of the electromagnetic field. Thus, the nonlinear term 
in (|3.1|) is not disturbing at this level, and the Poynting vector P preserves 
its meaning. 

By noting that J • J = 1 and that E • J = 0, we get another interesting 
relation: 

ar = a(E_j) _ m = _ c2(curlB) .j + cdivE (36) 

dt dt at 

Finally, one has: 

d 2 B , <9E 



curl—— = — c curl(curlB) + c curl(pJ) 



dt 2 dt 

= - c 2 V(divB) + c 2 AB + c curl(pJ) (3.7) 

from which we deduce the following second-order vector equation with a 
nonlinear forcing term: 

<9 2 B 

— - = c 2 AB + c curl(pJ) (3.8) 
ut 

that generalizes the second equation in (|2.20|) . We are sorry to announce 
that the "wave" equations for the fields E and B are no longer true. On the 
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other hand, it has emerged in section 2 that, in vector form, they are only 
a source of a lot of trouble. 

In the classical Maxwell equations the role of the field E can be inter- 
changed with that of field cB. This is not true for the new formulation. 
We will later see, in section 9, how to solve this problem. For the moment, 
we keep working with (|3.1|) - (|3.2j) - (|3.3|) . just because the theory will be more 
easy. In the coming sections 4 and 5, we will see how elegantly it is possible 
to solve the problems raised in section 2. 



4 Perfect spherical waves 

In the case of a plane wave of infinite extension, for both the Maxwell model 
and the new one, we are able to enforce the condition divE = and realize 
the orthogonality of the Poynting vectors with respect to the propagation 
fronts. Concerning a "fragment" of plane wave, the classical method runs 
into problems. However, with the new approach the situation radically 
improves. Let us see why. 

With the same assumptions of section 2, let be the square [0, 1] x 
[0,1]. We already noted that, on the sides {0}x]0, 1[ and {l}x]0, 1[, the 
quantities divE and curlB become infinite. Nevertheless, when substituted 
into equation 1)3.1)) . they come to a difference of the type +oo — oo. The two 
singular terms reciprocally cancel out, leaving a finite quantity, so that the 
equation has a chance to be satisfied. To show this, we can create a layer 
around the boundary of f2. Then, without developing singularities, we pass 
to the limit for the width of the layer tending to zero. The trick now works, 
because, in contrast to the classical Maxwell case, equation (|3.1|) can be 
satisfied exactly in all the points, since it is compatible with the condition 
divE 7^ 0. In the limit process we can also guarantee that the Poynting 
vectors remain parallel to the z-axis. Therefore, the fragment does not 
change its transversal shape. Explict computations will be carried out in 
section 5, in the case in which Q is a circle. 

For example, the situation represented in figure 2 is perfectly allowed for 
by our equations, except near the lower and the upper sides. Actually, on 
the sides ]0, l[x{0} and ]0, l[x{l}, given that divB and curlE are singular, 
we still have problems (clearly equation 1)3.2)1 and 1)3 . 3 j) are not true) . Similar 
problems are encountered by modifying the polarization of the wave. These 
troubles will be solved in section 9, by unifying 1)3.2)1 and 1)3.3)1 in a single 
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equation similar to Q3.1JI . in such a way that the roles of E and cB are 
interchangeable . 

The case of a spherical wave is very interesting. Let us consider the 
transformation of coordinates given in (|2.6j) . Let us also suppose that the 
fields are given as in (|2.7[) . with u, v, w not depending on 6. We have: 

E x B = (uw, - uv, 0) (4.1) 



E x B s(u) , s , 

? (w, -v, 0) (4.2) 



|ExB| Vu 2 + w 2 
where s(u) = u/\u\ is the sign of u. 
The new equations in spherical coordinates become: 

* --(*+") + ? < 4 - 3 > 



c 7 
7l l 



COS ( 

u— — 

r v sm< 



+ Uc 



c s(u) 



2 1 cos 4> 

v r + -v + -ws H : — -to 

r r rsinr 



Vw 2 + 



10 



(4.4) 



2 1 cos ^ 

TV + -V + -«;«* H ; -W 

w t = -c 1 - + u r ) + cs(u) ^— v (4.5) 

V r > yjv 2 + w 2 

These expressions may seem rather complicated, but there is nothing to be 
afraid of. 

To avoid discontinuities, we also impose the boundary conditions (|2,14|) 
and (|2.15jl . Now, by choosing v = and w = cu, one obtains: 

ExB , „ /cu cos d>c \ , , , 

— — -divE = ( _^ + -^, 0, 4.6 

|E x B| \ rsmcp r r / 

Therefore, one gets: 



u t 



- K + 7) ^ 

v t = (4.8) 
(«r + z) (4-9) 



^ = - c 

r 
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Once again, the first and the last equations are equivalent, providing the 
general solution ()2.16|) . Anyway, this time, thanks to the nonlinear corrective 
term of the second equation is compatible with v = 0. We are not 

obliged to choose /(</>) = 1/ sin 0, in order to enforce the condition divE = 0, 
because this constraint is no longer required. The conclusion is that perfect 
spherical waves are admissible with the new model. The functions / e g 
may be truly arbitrary (the only restriction is /(0) = /(vr) = 0). Continuing 
with our analysis, we will construct later infinite other solutions which are 
unobtainable with the classical Maxwell model. 

In the perfect spherical case, the Poynting vector P = (cu 2 , 0, 0) only 
has the radial component different from zero. As expected, this component 
has a constant sign (even if u and w oscillate). Since the set of equations 
is of a hyperbolic type, we can introduce the characteristic curves. In the 
example of the spherical wave, such curves are semi straight-lines emanating 
from the point r = 0, and the vector J = P/|P| = (1, 0, 0) is aligned with 
them. The nonlinear term introduced in (j3.1j) does not adversely affect the 
behavior of the wave, because, with v = and w = cu, the corresponding 
equations (|4.7|) and Q4.9JI are linear. Therefore, the superposition principle is 
still valid. Any piece of information, present at the boundary of the sphere 
of radius r > 0, propagates radially at the speed of light, without being 
disturbed (except by the natural decay in intensity). The nonlinear effects 
of the model are latent. They show up when we try to force, with some 
external solicitations, the Poynting vector not to follow the characteristic 
lines. This circumstance will be taken into account in sections 7 and 8. 

In a very mild form, we can state that the divergence vanishes, by ob- 
serving that, for any T: 



that is divE has zero average when integrated over a period of time. Never- 
theless, in section 13, we will get an astonishing result. We will see that an 
electromagnetic wave produces, during its passage, a modification of space- 
time. In the new geometry, the 4-divergence of the electric field is zero. 
This could make it difficult, or even impossible, to set up an experiment 
that emphasizes the condition divE 7^ at some point. The measure could 
be affected by the modified space-time geometry in such a way the condition 
cannot be revealed. 

Among the stationary solutions we find: 




divE dt = 



(4.10) 



u(t,r, 4>) 



w(t,r, 4>) 



K 3 



(4.11) 



r sin 



r sin 
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as well as: 



u(t,r, cf>) 



v(t, r, (f)) = K<i cos 4> w(t, r, 



4>) 



K 2 s\n<j) (4.12) 



r sm 



with K\, K2, K3 arbitrary constants. In particular, we recognize the classical 
stationary electric field: 



whose divergence is zero for any r > 0. Unfortunately, most of these solu- 
tions show singularities. 

Due to the nonlinearity of the equations, the study of the interference of 
waves looks quite complicated. As long as the waves are such that divB = 
and divE = (as in the plane case) there are no problems, since the 
nonlinear terms do not actually activate. For waves of different shape, the 
situation may be truly intricate. In first approximation, the nonlinear effects 
should attenuate faster than the amplitude of the waves. Thus, at a certain 
distance, these anomalies may not normally be observed. Although we do 
not wish to discuss it here, the subject is of crucial relevance and deserves 
to be studied in more detail. 



5 Travelling signal-packets 

In this section, it is convenient for us to express our new set of equations 
in cylindrical coordinates. After taking (x,y,z) = (r cos 0,r cos 9, z), we 
assume that the fields are of the form B = (0, 0, u), E = (v, w, 0), where the 
first component is referred to the variable r, the second one to the variable 
z and the third one to the variable 9. Moreover, for simplicity, the functions 
u, v and w will not depend on 9. In cylindrical coordinates, the counterparts 
of equations (|Q |) -(|4^i j) -(|43 |) are: 



E = (K 2 r-\ 0, 0) 



u t = -w r + v z 



(5.1) 



V 

iV H 1- Wz 




(5.2) 




v r -\ h w z 



(5.3) 
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By setting v = 0, we can easily find solutions when u and w do not depend 
on z. In this case one has divB = and divE = 0. From l)5.1j) and ()5.3|) it 
is easy to get the equations: 



U tt = <? ^ + U r J W t t = <? + W r j (5.4) 

whose solutions are related to Bessel functions. 

Anyhow, extremely interesting solutions in cylindrical coordinates turn 
out to be the following ones: 

u(t, r, z) = g(t ± z/c)f(r) v(t, r, z) = ±c g(t ± z/c)f(r) 

w(t,r,z)=0 (5.5) 

Note that the divergence of E is equal to p = v r + r~ 1 v, so that it is 
different from zero, unless / is proportional to 1/r. The functions / and g 
can be arbitrary (to guarantee the continuity of the vector fields, we only 
impose /(0) = 0). The relations in (|5.5[) give raise to electromagnetic waves 
shifting at the speed of light along the z-axis. If / and g vanish outside a 
finite measure interval, for any fixed time t the wave is constrained inside a 
bounded cylinder. The packet travels unperturbed for an indefinite amount 
of time. The corresponding field E is perfecly radial and the vector J is 
parallel to the z-axis. 

Given ro > 0, suppose that / is zero for r > r$. Suppose also that /, 
in a neighborhood of ro has a sharp gradient. It is evident that the vector 
c 2 curlB — c(divE)J = (c 2 u^,0,0) = (±cg'f, 0, 0) remains bounded even if 
we let the derivative of / tend to oo at ro. Therefore, as we anticipated at 
the beginning of section 4, we can give a meaning to equation (|3.1|) . even if 
/ is a discontinuous function in ro- 

We can get a transport equation for the unknown p = divE by using the 
equation (|3.4|) . i.e.: 

dp 

— = — c div(pJ) = — cp divJ — cVp • J (5-6) 

which, thanks to the fact that J is a constant field, takes the simplified form 
of: 

at = ±c d~ z (5 - 7) 

with the sign depending on the orientation of J. 
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We recall that the Maxwell equations do not allow for the existence of 
solitary waves, as the ones we have introduced right now. Therefore, here 
we obtained another important result. 

The energy £ of these solitary waves is obtained by integrating the energy 
density, given by: ^eo(|E| 2 + c 2 |B| 2 ). Thus, one gets: 

r+oo r+oo 

£ = 2vre c 2 / f{r)rdr / g*(£)d£ (5.8) 

JO J-oo 

Suppose that, at an initial time to, the electromagnetic fields are as- 
signed compatibly with (|5.5|) . The vector J turns out to be automatically 
determined, then the wave is forced to move in the direction of J at speed 
c. There are no stationary solutions, unless g is constant. But, in this last 
case, due to (|5.8(l . the energy is not going to be finite. The wave-packets 
take their energy far away, with no dissipation, until they react with other 
waves or more complicated structures (like, for instance, particles). 

Let us study more closely the expressions given in (|5.5jl . Assume that, 
for r > 0, the function / is non negative, and that /(0) = 0. If the function g 
has a constant sign, we distinguish between two cases, depending if the sign 
is positive or negative (see figure 3). The sign determines the "orientation" 
of the vector curlB (note however that curlB is not exactly parallel to the z- 
axis, despite what is shown in figure 3). Then, we have subcases, depending 
whether E is directed toward the z-axis or not. In conclusion, two possible 
types of solitary waves can occur, depending on the orientation of the electric 
field (external or internal). These will be denoted by 7 + and j~ , respectively. 
In figure 3, J indicates the direction of motion. Of course, g could also have 
a non-constant sign. In this case, the corresponding wave can be seen as a 
sequence of waves of type j + and j~ , shifting one after the other. 

In vacuum, the electromagnetic fields at rest, are assumed to be identi- 
cally zero. During the passage of a soliton, the calm is momentarily broken 
only at the points "touched" by the wave. The information shifts, but does 
not irradiate. Thus, also if the wave-packet displays a negative or positive 
sign, depending on the orientation of E, this is not in relation to what is 
usually called electric charge. Hence, as long as the cylinders containing two 
different solitons do not collide, they can cross very near without influencing 
each other. On the contrary, we expect some scattering phenomena, through 
a mechanism that will be studied in section 8. 

If we place a mirror parallel to the z-axis, at some distance from it, 
the reflected image of the travelling wave-packet will be the same wave- 
packet shifting in the opposite direction (because curlB changes sign, while 
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Figure 3: Behavior of the electric held for the two possible wave-packets (section 
for a fixed angle 9): 7~ shifting from left to right, j + shifting from right to left, 
7~ shifting from right to left, 7 + shifting from left to right. 

E maintains its orientation). This disagrees with our common sense. In 
other words, equation (|3.1|) does not preserve plane symmetries. The same 
is true for the Maxwell equations. In both cases we have no elements to 
decide the correct sign of the vector product x (left-hand or right-hand). 
As a matter of fact, without modifying the equations, a change of the sign 
of x can be compensated for by a change of the sign of the electric (or the 
magnetic) field. To learn more about this problem we need to study the 
interactions between waves and matter. Hence, for the moment, we have 
no sufficient information to understand from which part of the mirror is our 
universe. An answer to these crucial questions will be given in section 8. 

In cylindrical coordinates, we can find many other solitonic solutions. 
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Here is another example: 



E = (±cu, 0, =fcv) B = (v, 0, u) (5.9) 

with u = fi(r, 9)g(t ± z/c) and v = fa(j, 9)g(t db z/c). In order to fulfill the 
condition divB = 0, it is necessary to impose: 

dfx d ( r f2) 

J± + S£l =a (5.10) 

Note that E • B = and J = (0, =Fl, 0). Except for the condition (f5~TU|) . the 
functions /i and /2 are arbitrary, so that the new solutions are very general. 
Actually, they include the previous ones (take fi = f and fa = 0). In section 
9, we will remove the condition divB = 0, allowing for the existence of even 
more solutions. We can force f\ and fa to vanish outside a bidimensional 
domain VL. Again, this plane front, modulated by the function g, travels 
along the z-axis at the speed of light. 

We may reasonably expect that these solitary solutions are modified, or 
even destroyed, when they encounter another external electromagnetic field. 
In fact, the equations being nonlinear, the superposition principle does not 
hold, in particular if the motion is disturbed in a way that is in contrast to 
the natural evolution along the characteristic curves. As far as we know, 
there are no documented experiments evidencing these facts. In section 
11, we instead examine the behavior of solitary waves under the action of 
gravitational fields. 

An electromagnetic radiation can be suitably considered as an envelope 
of solitary waves, travelling in the same direction. If, transversally, these 
solitons are of infinitesimal size, they can be compared to "light rays" . This 
observation clarifies how a wave can be viewed at the same time as a whole 
electromagnetic phenomenon and as a bundle of infinite microscopic rays. 

And then there are photons. They are also pure electromagnetic manifes- 
tations, but, unfortunately, they are not modelled by the Maxwell equations. 
Modern atomic and subatomic physics would not exist without photons, yet 
they find no place in the classical theory of electromagnetism. This is an 
unpleasant gap. Although physicists are acquainted with this dualism, the 
general framework remains blurred. From our new standpoint, we contend 
that the photons observed in nature have very good chances to be modelled 
by the equations introduced here. As a matter of fact, we have enough free- 
dom to be able to build solutions (no matter how complicated) resembling 
real photons. We can assign a frequency to them, longitudinally or transver- 
sally. Then, we know that they must move at the speed of light and can 
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have finite energy, given by the energy of their electromagnetic vector fields. 
They can be "positive" or "negative" without being electrical charges. Even 
with no mass at rest, their motion can be distorted by gravitational fields 
(see section 11). A concept of spin can be also introduced (see section 15). 
If what we are proposing here is a new functioning model for electromag- 
netism (and we will collect other evidence supporting this hypothesis), then 
it explains why photons can be self-contained elementary entities and elec- 
tromagnetic emissions at the same time. In this case, a first link between 
classical and quantum physics is set forth. 



6 Lagrangian formulation 

In order to recover the equations (|3.1|) - ()3.2|) - (|3.3j) from the principle of least 
action, we follow the same path bringing to the classical Maxwell equations. 
Thus, we introduce the scalar potential <& and the vector potential A = 
(A\, A2, A 3 ), such that: 

1 1 dA 
B = - curlA E = V$ (6.1) 

c c dt v ; 

From the above definitions we easily get the equations: divB = and 
J^-B = — curlE. The third equation (|3.1(l is going to be deduced from the 
minimization of a suitable action function. 

Let us first note that, by taking the pontential equal to zero and 
setting £ = t — (xJ\ + yJ% + zJ 3 )/c, one obtains: 

1 dA OA 1 dA 

E = --— cB = V£x— = -Jx-- 6.2 

c ok, ot, c a£ 

This allows us to infer that B is orthogonal to E and that |E| = c|B| (see 
also QUI, p. 126). 

Successively, for i and k between and 3, we introduce the electromag- 
netic tensor: 

:) 1, n \ 

(6.3) 





dA k 


dA, 


Fik — 


dxi 


dx k 


X2,X 3 ) 


= {ct, 


-x,-y, 
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-Ex 


—E2 


Ei 





-cB 3 


E2 


cB 3 





\E 3 


-cB 2 





where Aq = $ and (xo, aci, x 2 , x 3 ) = (ct, —x, —y, —z). Explicitly, we have: 
E.i. - 7' r, \r (6-4) 



-E 3 \ 
cB 2 
-cB x 
/ 
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with E = (Ei,E 2 ,Ez) and B = (B 1 ,B 2 ,B 3 ). Replacing E by -E, one gets 
instead the contravariant tensor F tk : 



( 
-Ex 

-E 2 

\-E 3 



Ex 


cB 3 
-cB 2 



E 2 E3 \ 

-cB 3 cB 2 

-cB L 
cB x / 



(6.5) 



Therefore, up to multiplicative constants, the action turns out to be (see for 
instance [S], p. 596): 



S 



F ik F ik dxodxidx 2 dx 3 = c I F ik F llt dxdydzdt (6.6) 



pik , 



pik 



where, summing-up over repeated indices, the Lagrangian is L = F^F 1 
2(c 2 |B| 2 — |E| 2 ). As customary, the variations are functions SAi, having 
compact support both in space and time (between two fixed instants). 
With well-known results, one obtains: 



5S 



dF 



ik 



dxi 



SA{ dxodxidx 2 dxz 



(6.7) 



Imposing SS = 0, the corresponding Euler equations are exactly the stan- 
dard Maxwell equations. As a matter of fact, due to the arbitrariness of 
the variations 5Ai, one recovers: — ^-E ik = (for i = 0,1,2,3), that is 

equivalent to write divE = and t|E = c 2 curlB. 

Let us now introduce a novelty. We require that the variations 5Ai are 
subjected to a certain constraint, so that the conclusions are going to be 
different. Actually, we impose the condition: 



5<S> 



J-5A 



5A - Jt 5Ai - J?5A 9 . - JzSA? 



'i< 







(6.8) 



where J is the vector (E x B)/|E x B|, already defined in section 3. The 
relation ()6.8|) says for instance that, if the vector variation (5A\, 5A 2 , 5A 3 ) 
locally belongs to the tangent plane generated by E and B, then the vari- 
ation 5Aq is zero. Although for the moment we can only provide vague 
explanations, we assert that ()6.8|) is the germ of the Huygens principle. The 
picture will become more focused as we proceed with our analysis. 

Consequently, we discover that the 4- vector —jj^-F %k is not identically 
vanishing, but, due to (|6.8j) . it must have a component along the 4- vector 
(1, — J). This leads to: 

divE = A (6.9) 
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1 <9E 

ccurlB = - AJ (6.10) 

c at 

where the parameter A is a Lagrange multiplier. By eliminating A we easily 
arrive at equation (|3.1|) . Thanks to (|6.8|) . the set of possible variations is 
smaller than the set in which we impose no restrictions at all. Therefore, the 
equation 5S = is now less stringent. As we already know, this shows that 
(|3.1|) admits a space of solutions which is larger than the one corresponding 
to (j2~Tj) and (j2~2]) together. 

Using the electromagnetic tensor, the equation (|3.1|) can be written as: 

( Qpik QpQk \ 

- c - — + — — Ji = fori = 1,2, 3 (6.11) 

By defining Jo = —1, the above relation is also trivially satisfied for i = 0. 
In section 11, within the framework of general relativity, we will be able to 
write (|6.11|) in invariant tensor form. 

Let us try to understand the reason for the constraint (|6.8|) . As far as 
the evolution of a plane or a spherical wave is concerned (and, surely, in 
more general cases), it is easy to check that: 

A = $J (6.12) 

A straightforward way to get (|6.12j) is from explicit solutions. For example, 
in spherical coordinates, we can use ()2.16j) in order to find: 



E 



(0, l/^t - r/c), 0) B = (0, 0, l/(^(t- r /c)) 



J = (1,0,0) A = (-F(<j>)g(t-r/c), 0, 0) $ = -F(<j))g(t-r fc) (6.13) 

where F is such that F' = f. The relation (|6.12j) is also true in the case of 
solitons. In fact, in cylindrical coordinates, thanks to IJ5.5|) one has: 



E = (±cf(r)g(t±z/c), 0, Oj B = (0, 0, f(r)g(t ± z/c) 
J = (0,=Fl,0) A = (o, cF(r)g(t±z/c), 6) $ = ^cF{r)g(t±z/c) (6.14) 



where F is such that F' = f and /(0) = 0. Note that, in general, $ and A 
are not uniquely determined. However, there exists at least one choice of 
and A such that (J6.12|) is satisfied. 
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Now, the equation A = <1?J implies |<J>| = |A|, or equivalently: <5 2 — 
| A| 2 = 0. Taking the variation of the last relation brings us to the constraint: 

2($ 5$ - A-<$a) =2$ (S<$> - J-8A\ =0 (6.15) 

which is the same as in (j6.8|) . Another way to get ()6.8|) is by directly evalu- 
ating the variation of (|6.12j) : 

5A = <5($J) = 5$ J + $ 5J J • SA = 5<S> (6.16) 

where we used that |J| = 1 and 53 ■ J = (a normalized vector is orthogonal 
to its variation). 

Obviously, the vector relation A = <J>J implies the scalar relation: 

$ = J- A (6.17) 

obtainable after scalar multiplication by J and by observing that |J| = 
1. Another way of expressing (J6.17|) is to require that the scalar product 
between the 4-vectors (3>,A) and (1,— J) is zero. This makes (|6,17[) an 
invariant in the context of general relativity. 

In conclusion, the equation (|3.1jl can be recovered from the constrained 
minimization of the action function associated with the classical Lagrangian. 
The constraint originates from (|6.12j) which says, in particular, that A is 
lined up with J. As will be better explained in section 10, such a condition 
characterizes the propagation of waves, whose evolution is ruled by the Huy- 
gens principle. From now on, these will be called "free waves". In sections 
7 and 8, we will see that not all waves are of this type. 

An interesting relation, that is directly obtained by checking (|6.13|) or 
()6.14|) . is the following one: 

E + cJxB = (6.18) 

The above equation is extremely important, since it represents another char- 
acterization of free waves, perhaps simpler than (|6,12|) , Indeed, it is the 
analogous of the Lorentz law for moving electric charges (recall that cJ is a 
velocity), even if here there are no particles around. As it will be explained 
in the coming sections, equation (J6.18|) tells us that the mechanical forces 
acting on a free wave are zero. Therefore, the wave actually moves without 
external disturbances and in agreement with the Huygens principle. We will 
prove all these statements in sections 9 and 10. 
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The covariant version of (|6.18|) is: 



- cB + J x E = (6.19) 

obtained by vector multiplication of (|6.18|) by J. Both ()6.18|> and ()6.19j) can 
be trivially deduced from the orthogonality of E with respect to B, and by 
the equality |E| = c|B| (see also the beginning of this section). Therefore, 
they have quite a general validity. 

Before ending this section, we will collect a few more properties. We 
begin by considering the following writing: 

|E| 2 -c 2 |B| 2 = E-(--^ - V<f>) - cB-curlA 
11 11 \ c dt J 

Id, . 1 8E 

— c A • curlB — c div(A x B) 

- div($E) + $ divE - c div(A x B) (6.20) 

where we used the definitions in (|6.1|) and known formulas of vector calculus. 
Introducing the constraint A = <3?J (hence also $ = J ■ A) we have that A 
is orthogonal to E, because so it is J. Thus, (|6.2U|) can be simplified and 
becomes: 

|E| 2 -c 2 |B| 2 = -div($(E + cJxB) 



1 / f)V \ 
+ -f — - c 2 curlB + c(divE)JJ - A (6.21) 

Then, it is interesting to point out that, when both (|6.18j) and (|3.1|) are 
satisfied, the Lagrangian vanishes. On the contrary, when E is orthogonal 
to B and |E| = c|B| (so that (|6.18|) holds), the relation (|6.21|) reveals that 
imposing equation ()3.1|) is a natural requisite. 

Finally, due to ()6.1|) . equation (|3. 1|) entails: 

1 <9 2 A ^9$ , A , (\ dA . \ , 

pt-tt - V— = c curl(curlA) + c - div^— + A$ J 

c ot at \c at J 

= -c AA + cV(divA) + ( + c A$\ J (6.22) 



28 



We also assume the following Lorentz condition: 



id® 

divA + -— = 

c at 



(6.23) 



which is known to be an invariant expression in general relativity. With the 
help of fl6~23]) . from (l?H2l we get: 



which is in perfect agreement with the continuity equation (|3.4[) . once we 
define / o = c _2 Jp-<I ) — A<£. Note however that (|3,4|) is true independently of 
(|6.23f) and (|6.24l) . In the classical Maxwell case, the right and the left terms 
of (|6.24|) both vanish, providing, together with (|2.20|) . additional "wave" 
equations. This is not necessarily true in our case. 



7 Encounter of a wave with an obstacle 

In this section and in the following one, we study what happens to a free 
wave when it meets an obstacle that we define of "mechanical type". The 
phenomenon can be extremely complicated, therefore the analysis will be 
carried out on very simple cases. For the moment, here we just discuss some 
basic facts, trying to catch the underlying ideas. In section 9, we formalize 
the problem better, by writing down the final equations. 

We first take into account an example concerning the reflection of an 
electromagnetic radiation. We assume to have a monocromatic plane wave, 
linearly polarized, which is totally reflected by a perfectly-conducting metal- 
lic wall. Hence, we suppose that there is no refraction at all. In Cartesian 
coordinates, the wall corresponds to the plane y = 0. Initially, the wave- 
front propagates forming an angle C 7^ with the y-axis. Referring to figure 
4, the incident wave is described by the fields: 




(6.24) 





(7.1) 



The reflected electric field E( r ) is such that, for y = 0, the component 
parallel to the obstacle, of the vector E^ r ^ + EW, vanishes (see for instance 
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Figure 4: Reflection of a plane wave when the magnetic field is normal to the plane 
of incidence x = 0. The vectors and are therefore orthogonal to the page. 

[2], p. 270). Concerning the magnetic field at y = 0, we have: E$( r ) = B®. 
Then, one easily gets: 

= ^0, — csin£, — ccosC) smuj(t + (ycos( — zsm()/c) 

B (r) = (l, 0, 0) sin w (t + (y cos C-£ sin C)/c) (7.2) 

To justify the imposition of the boundary conditions on y = 0, it is cus- 
tomary to assume the existence of instantaneous electric currents on the 
conducting surface, which force the tangential electric field to zero (see [3], 
p. 558, and [§], p. 335). 

After reflection, the wave is very similar to the incident one, with the 
exception that = (0, cos £, sin Q has changed in j( r ) = (0, — cos C, sin Q. 
Note that the relation (|6.18|) is valid both for the incident and the reflected 
waves. Nevertheless, during the impact, in which an instantaneous flip of 
the signs occurs, the wave-front at y = does not show a behavior con- 
sistent with the one corresponding to a free wave. Each wave-front, when 
reaching the wall, evolves for a single moment without respecting the eikonal 
equation. Of course, what we are considering here is just a mathematical 
idealization. More realistically, the wall is made of matter, so it would be 
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Figure 5: Reflection of a plane wave when the electric field is normal to the plane 
of incidence x = 0. The vectors EW and E( r ) are therefore orthogonal to the page. 



more correct to check what happens to the wave when it interacts with the 
atoms of the wall. Anyway, we do not think it is useful to study this more 
complicated situation, since it only modifies the form but not the substance 
of facts. We believe that the main idea has already been outlined: when 
encountering an obstacle a free wave may lose its likeness and become, for 
a small amount of time, a "constrained wave" . The reaction of the obstacle 
can be so strong that, as in the case of the reflecting wall, the wave- fronts 
are forced to modify abruptly the direction of their movement. What we 
would like to do in the following pages, is to understand what happens at 
those instants. 

We recall that, for y = 0, the vector E( r ) + does not have the same 
length of EW before the impact (or E( r ) after the impact). Therefore, some 
electromagnetic energy turns out to be missing. We conjecture that it has 
been spent, with the help of the mechanical constraint, to allow the variation 
of the direction of the wave- front of an angle n — 2(. For a moment, this 
energy has gone elsewhere, compensated by something which is not of an 
electromagnetic kind. We would like to find out what is. To this end, let us 
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introduce a new vector field: 




if y ^0 



G 



j(r) _ 



(7.3) 



if y = 



Note that, dimensionally, the vector G corresponds to an acceleration. In 
the case under study, at y = 0, G is oriented as the vector (0, —1,0), with 
an infinite magnitude. If we imagine the wave as a bundle of rays (we saw in 
section 5 that the two things are strictly related) , then G provides a measure 
of their curvature. When we are dealing with a free wave, the rays proceed 
along straight-lines. Corresponding to this situation, we have G = 0. When 
the rays hit the wall, then G becomes different from zero. In the particular 
case we are examining, G is a singularly distributed field, but if we allow the 
rays to change their trajectories smoothly, then G is going to be finite. This 
would give the idea of a centripetal time-varying acceleration, responsable 
for the rotation of the rays and the corresponding wave-fronts. 

In the coming sections, we will make evident that a nonvanishing vector 
G appears each time a wave evolves without following the Huygens principle, 
as a consequence of external perturbations. Some suitable way of estimating 
the magnitude G should allow us to compensate the missing electromagnetic 
energy. By the way, the theory is not going to be easy. In section 13, we will 
discover that, in order to change the trajectories of the rays, it is necessary to 
pass through a modification of the space-time structure. Thus, we cannot be 
more precise, until we are ready to carry out our analysis in the framework 
of general relativity. Before that, we have to work a little more on the 
definition (|7.3|) . This will be done in section 9. For the moment, let us put 
together other basic facts. 

We now vary the polarization of the incident wave of 90 degrees (see 
figure 5). For y = 0, the component, orthogonal to the reflecting wall, of 
the magnetic field must be zero. At the same time, the whole electric field 
vanishes (since, for y approaching to zero, the field E( r ) is opposite to EW). 
Therefore, more electromagnetic energy is missing when the wave hits the 
wall. Actually, in this case, the effects of the obstacle are stronger: together 
with the deviation of the wave- front, we also observe a torsion that modifies 
the polarization by 180 degrees. Such a torsion process is instantaneuos, but 
qualitatively similar to that corresponding to a circularly-polarized plane 
wave, like the one for example expressed by the following fields: 



E = (ccosuj{t — y/c), 0, csintvit — y/c)) 
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B = (sinw(t - y/c), 0, - cos oj(t - y/c)) (7.4) 

where the polarization constantly changes at finite speed. We note that this 
wave is also a solution to Maxwell equations and is more "energetic" than 
the corresponding one with constant polarization. 

If, together with reflection, we also have refraction within a medium of 
different density, the study is more involved. When reaching the plane of 
reflection the rays bifurcate. In our opinion, this is due to the arbitrariness 
of the vector J at the time of the impact (for example because P is zero). 
The incident wave splits into two solutions (the reflected and the refracted 
waves), both compatible with the same boundary conditions on the plane 
y = 0, brought by the incoming solution. We do not necessarily have bifur- 
cation each time that P is zero (this actually happens very frequently), but 
only when, in addition to this, suitable uncommon conditions are verified. 

Based on the above observations, we conclude this section with a little 
philosophical dissertation. The equations (|3.1|) - (|3.2j) - ([3.3|) are of local and 
deterministic type. Hence, for given initial data, the solution is uniquely 
determined. Nevertheless, there may be circumstances in which, following 
the evolution of a certain solution, several other branches of solutions of the 
equation (|3.1j) may be admissible. As we noticed, this could be true because 
the evolution of the vector J turns out to be compatible with different op- 
tions. Thus, the original solution splits, and this event is also deterministic. 
As a matter of fact, when the right conditions are fulfilled, an incident wave 
bifurcates, giving raise to a reflected and a refracted wave. There is no 
uncertainty: both the solutions are systematically observed. However, this 
situation becomes extremely unstable when reversing time. We are indeed 
allowed to think that, marching backward in time, the inverse of a bifurca- 
tion phenomenon could occur. In this case, as a film runs in reverse, two 
waves would meet in perfect phase and melt, producing a single wave. Nev- 
ertheless, this event has zero probability of happening. A little perturbation 
is sufficient to modify completely the evolution of the two waves, with no 
chances of seeing their fusion. 

In conclusion, our equations are of hyperbolic type, deterministic and en- 
ergy preserving. Nevertheless the particular nature of the nonlinear term can 
be a source of instabilities when reversing the sign of time. The consequence 
is that some original information can get lost, and there is no practical way 
to recover it by following the reverse path. We ask ourselves if this could be 
an explanation (at least in part) for the second law of thermodynamics. 
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8 Diffraction phenomena 



We continue our qualitative analysis on the interaction of waves with mat- 
ter. The second example that we take into account is the developing of 
diffraction, where an electromagnetic wave encounters an aperture. Once 
again, for simplicity, the device used to generate the phenomenon is a pure 
mathematical abstraction. So, let us suppose that a plane monocromatic 
wave propagates in the direction of the y-axis (with y increasing) and hits 
a perfect-conducting wall at y = 0. This time, however, there is a passage 
through the strip < z < a, for some positive width a. 

For y < 0, the incident wave is polarized as in ([7.1JI with £ = 0. As a 
consequence, we have: 

E = (0, 0, csmuj(t-y/c)) B = (smuj(t - y/c), 0, 0) 

ExB = P = (0, c[smuj{t-y/c)] 2 , 0) 

J = — — = (0, 1, 0) divE = divB = (8.1) 

For z < and z > a the wave is reflected. Actually, it is not a perfect reflec- 
tion, since it is affected by some perturbations originating at the boundaries 
of the aperture. Here, we neglet this aspect and focus our attention on the 
study of the sources of these disturbances. 

We assume that along the two straight-lines y = 0, z = and y = 0, 
z = a, instantaneous electric currents push the electric field to zero, so that 
the condition E = is enforced (see (Hj, p. 559). At the instant in which 
the wave reaches the obstacle, the electric field develops discontinuities. 
Therefore, its divergence is a concentrated Dirac distribution along the above 
mentioned straight-lines. Thus, for y = 0, it is easy to realize that: 

divE = c [5 Q (z) - 5 a (z)]smujt (8.2) 

On the other hand, for < z < a, we can expect that: 

lim P = (0, c(sinwt) 2 , 0) (8.3) 

because these are points in which the wave does not hit the obstacle. 

Next, let us examine more in detail what happens on the straight-lines 
y = 0, z = and y = 0, z = a. Using equation (|3.1|) and neglecting the term 
cur IB, the remaining nonlinear term brings an instantaneous rotation of the 
electric field. As a matter of fact, let us define in Cartesian coordinates 



34 



E = (0,v,w) and B = (u, 0, 0), where u, v and w do not depend on x. 
When y < 0, we trivially have v = 0. In terms of the new unknowns, the 
equation Q3.1JI takes the form: 



the sign of 

We have u z = 0. Neglecting u y (which is bounded), the system (|8.4j) 
is equivalent to a rotation with angular speed H (note that, for z = and 
z = a, 3 is infinite because so it is i« 2 ). The rotation is anti-clockwise at the 
points (cc, 0,0). It is clockwise at the points (x,0, a). This is true whatever 
the sign of the electric field (note that E and B change sign together and 
J always maintains the same orientation). Even if the rotation is at infinite 
speed, the rotation angle is finite and depends on the magnitude of the 
fields. In practice, by equation (|3.1j) . a sudden change of the quantity divE 
is balanced by a modification of JrE, forcing the field to vary in the direction 



This behavior is not at all in agreement with what is usually observed. 
The diffraction is a diffusive phenomenon. Therefore, referring to figure 6, we 
should expect a clockwise rotation at (x, 0, 0) and an anti-clockwise rotation 
at (x,0,a). We can correct this inconsistency by changing the sign to the 
vector product x, which amounts, in other terms, to inverting the sign of the 
electric field. The conclusion we can draw is quite astonishing. According 
to our equations, the standard right-handed product x is not suitable for 
describing a natural event like diffraction. Answering the question raised in 
section 5, the correct side of the mirror is the one where x is left-handed. 
This means that, if we do not wish to modify the model equations, we need to 
switch the polarity of the electric fields, in such a way an electron turns out 
to be positive and a proton negative. We will return to this subject in section 
15. The asymmetry of our universe and, consequently, the determination of 
its parity, is a problem that emerged about 50 years ago. The effects of this 
dichotomy were predicted by Lee and Yang (see for instance ^21 > P-534), 
but the reasons for preferring left or right have still to be found. If the 
above arguments are free from errors, then here may lie the solution to the 
problem. 

In a more realistic situation, the change of direction of the Poynting 
vector field is not instantaneous, but develops smoothly. It is very important 
to point out that the diffraction comes as a consequence of the nonlinearity 
of the model equation (j3.1|h and not from the imposition of some artificial 
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Figure 6: Plane wave encountering a reflecting wall with an aperture. Qualitative 
behavior of the Gelds E and 3 just before the impact and during the impact. The 
vector B is orthogonal to the page. 

boundary conditions at the points (x, 0, 0) and (x, 0, a), as supposed by other 
theories (see for instance chapter XI). The successive evolution of the 
wave after the obstacle follows the Huygens principle. In fact, for y > 0, the 
wave is free. It displays a slight diffusion due to the rearrangement of the 
electromagnetic fields, described before. It is well-known that the behavior 
depends on the ratio between the width a and the frequency u/2it. We do 
not investigate this aspect, since it has been intensively studied in the past. 
Here, we were only concerned with detecting the mechanism that leads to 
the deflection of the rays, when they meet the border of the aperture. The 
same mechanism takes place in the scattering of two solitons (see section 5), 
when their electromagnetic fields collide, with reciprocal influence. 

Some quantitative information can be recovered by examining equation 
(|3.5j) . We know that the Poynting vector P presents a natural pulsating 
variation along the direction of propagation of the wave. But, during the 
encounter with the obstacle, we have to add another variation, due to the 
instantaneous transversal change to the flux of energy. As we can see from 
figure 6, during the impact, the vector field P shows some dispersion and its 



36 



divergence suddenly grows. Thus, for just a moment, the quantity — c 2 divP 
registers a negative peek. To restore the normal energy flux in ()3.5|) . some 
corrective term, taking care of the "reaction" of the obstacle, should be 
added. In fact, as in the previous section, the change of curvature of the 
rays is accompanied by the creation of a new vector field G, concentrated 
on the obstacle. We better formalize this idea in section 9. 

Anyhow, in spite of the good achievements, we still have some problems 
to fix. Suppose that the incident wave is polarized in a different way, for 
instance by applying a rotation of 90 degrees: 

E = (cBanu(t-y/c), 0,0) B = (0, 0, -emu(t-y/c)) (8.5) 

On the contact with the straight-lines y = 0, z = and y = 0, z = a, we 
should now observe a prompt change of the magnetic field B, along the 
direction z. This event is not modelled by our equations, since we need to 
suppose that divB can be different from zero. In order to proceed, it is 
necessary to correct the model in such a way that the fields E and cB have 
the same role, as in the classical Maxwell equations. We also discuss this in 
the coming section. 



9 Adding the mechanical terms 

Based on some problems emerged in the previous sections, we make a first 
adaptation of the set of equations (|3.1|) - (|3,2|) - H3.3|) . with the aim of bringing 
to the same level the two fields E and cB. Thus, the new formulation is: 

<9E „ E x B . . 

— = c 2 curlB - c divE — — — (9.1 
dt |E x B| v ' 

$B E x B 

— = - curlE - c divB r (9.2) 

dt |E x B| y ' 

Now, H9.1JI - H9.2JI do not change if we replace E by cB and cB by — E, as 
in the Maxwell equations. This makes the model invariant under change 
of polarization. Clearly, if the divergence of B vanishes, we come back to 
equations H3.2J) - H3.3|) . The possibility for divB to be different from zero, does 
not imply the existence of magnetic monopoles, as the condition divE ^ 
does not imply the existence of electrical charges. However, the issue is 
delicate, and will be further discussed at the end of section 14 and in section 
15. 
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The modified version Q9.1|) - I)9.2JI allows for an even greater space of so- 
lutions. The spherical waves analyzed in section 4 can be now constructed 
with the electric field following the parallels, and the magnetic field following 
the meridians. In this case, we have divE = 0. Different other intermediate 
polarizations, which may also vary in time, can be considered as well. 

Finally, we can also get solutions as the one shown in figure 2. It is 
sufficient to set: 

E = (cf{x,y)g(t-z/c), 0, 0) 

B = (0, f{x,y)g(t - z/c), 0) J = (0, 0, 1) (9.3) 

where / is an arbitrary function, decaying to zero near the boundary of a 
bidimensional domain £1. The function g can also be arbitrary. If g vanishes 
outside a finite interval, then the solution in (|9.3|) is a travelling soliton, as 
the ones considered in section 5. The difference with respect to section 5 is 
that those solutions satisfied the condition divB = 0. Hence, in cylindrical 
coordinates, the solutions given in ()5.9)) satisfy (|9.1j) - ()9.2|) . without the need 
to enforce relation 1)5.10(1 . 

However, we are not completely satisfied yet. There is too much symme- 
try now, while we know that, in most natural phenomena, the role of fields 
E and cB is very well differentiated. Actually, the difference is detectable 
when a wave interacts with matter. We can take for example the case of the 
wire-grid polarizers, where an incident wave hits a grate of parallel metallic 
wires. If the wave is polarized with the electric field orthogonal to the wires, 
then it passes the obstacle almost undisturbed (if its wave-length is much 
smaller than the distance between two wires of the grate). If the electric 
field has a component along the direction of the wire, then the wave changes 
the polarization of a certain angle. 

Insisting on a similar example, we can go back to the end of section 
8. We are now able to study the diffraction of the wave given in 1)8.5(1 . 
where the electric field is parallel to the x-axis. Using the equations 1)9.1(1 - 
1)9.2(1 . we come to the same conclusions obtained for the wave given in 1)8.1(1 . 
polarized in another way. But this result is incorrect, because in the case of 
the wave 1)8.5)1 . together with the diffraction of the rays, there should be a 
change of the polarization after passing the obstacle, which is not modelled 
by the equations ()9. 1(1 - 1)9.2(1 . and which is not present in the case of the wave 
polarized as in (|8.1|) . 

In addition to the above observations, we also note in the reflection- 
refraction phenomenon, that the way the incident wave is polarized affects 
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the final result. Thus, it is necessary to further improve the model. For this 
purpose it will be useful the material collected in sections 6, 7 and 8. 

We need to introduce new vector fields (not of electromagnetic type), 
which are activated each time a free wave becomes a constrained wave. Let 
us begin by defining a velocity vector field V. We will ask all the vectors to 
be of constant norm, in particular: |V| = c, where c is the speed of light. 
Therefore, what really matters is the orientation of the vectors. The idea 
is that V is the tangent vector field to a bundle of light rays. An example 
is given by the vector field cJ, introduced in section 3, representing the 
direction of propagation of a wave- front. 

Afterwords, we propose the following system of time-dependent partial 
differential equations, with three unknown vector fields: 



<9E 
~dt 



(fcurlB - (divE)V (9.4) 



<9B 

— = -curlE - (divB)V (9.5) 

dV 

= _(V-V)V + |i(E + VxB) (9.6) 

at 

The constant fi is dimensionally equivalent to an electric charge divided 
by a mass. Finally, we add the condition previously anticipated: 

|V| = c (9.7) 

Concerning the choice of the norm in (|9.7|) . the discussion is postponed to 
section 15. 

It is customary, in fluid mechanics, to introduce the material (or sub- 
stantial) derivative: 

DV dV 

G ~ ~Dt ~ ~dt + (V ' V > V (9 ' 8) 

where G is an acceleration field. Hence, the equation (|9.6j) is equivalently 
written as: 

DV 

— = M (E + V x B) (9.9) 

Geometrically, the vector -j^V provides a measure of the curvature of the 
stream-lines, which in this case are identified with the rays (recall (|7.3j0 . As 
will become clearer, the knowledge of the vector field V is secondary with 
respect to the determination of its variation G. 
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Again, we assume to be in vacuum, with no particles of any kind around. 
In spite of that, the equation Q9.9JI has strong similarity with the Lorentz 
law for a density of charge moving at the speed of light. Actually, all the 
ingredients are there. Multiplying by a mass, the left-hand side in Q9.9[) 
is a force: its component along the direction of motion turns out to be 
proportional to the electric force field, while the orthogonal component is 
proportional to the magnetic force field. As we can see, the symmetry is 
broken, so that the electric and the magnetic fields cannot be interchanged 
anymore. But, this only happens in the case of constrained waves (G 7^ 0). 
For free waves, we recall that the relation (|6.18|) . corresponding to V = cJ 
and G = 0, is true. Moreover, replacing E by cB and cB by — E, we obtain 
the relation ()6.19j) . which is also true. Therefore, all the free waves, no 
matter what kind of polarization they have, are included in the new model. 
The interesting part is to study the behavior of constrained waves. We will 
discuss some general properties in section 10. 

Before going ahead, we feel that some clarification is necessary concern- 
ing the meaning of the word "mass" , used, perhaps improperly, several times 
in the paper. In our discussion, there are no masses in classical sense, since 
there are no elementary particles of any sort. Nevertheless, we needed to 
make distinction, in terms of dimensionality, between electromagnetic and 
mechanical (later they will be called gravitational) fields. This responsabil- 
ity has been given to the constant //, which provides the dimensional link 
between the two "flavors". Although other names could have been appro- 
priate to this purpose, the choice of the term "mass" is not incidental, since, 
as we proceed with our arguments, it will come out to be consistent with 
the standard setting. 



10 Properties of the new set of equations 

The new system of equations (|9.4)) - (|9.5[) - (|9.6)) is able to describe electro- 
magnetic phenomena where the wave- front, locally evolving in the direction 
determined by V, could be subjected to transversal perturbations modifying 
the trajectories of the rays. The propagation of the wave is governed by the 
first two equations. Through a feed-back process, the third equation, from 
the current knowledge of the local electromagnetic fields, allows for the de- 
termination of 157 V, setting up the new direction of motion. This coupling 
is possible because we have been able to include the vector V in the descrip- 
tion of the electromagnetic part, in the same way the term cJ was added to 
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(|3.1I) , Thus, we got a remarkable result: a link between electromagnetic and 
mechanical forces. Using the standard Maxwell equations such a connection 
could never be established. 

Let us continue with our analysis. From known formulas of vector cal- 
culus, we first deduce that: 

/,V . V = ^ ■ V - (-VIVI 2 - VxcurlvVv = (10.1) 

dt V2 1 1 / v ' 



Dt dt 

7lA7"|2 n nnrl 

at 



where we used that V|V| 2 = and that J^V is orthogonal to V, since V 
has constant norm. 

We recall that, by definition, E- J = 0. Similarly, by (|1().1|) and by scalar 
multiplication of (|9.9|) by V, one easily gets: 

E-V = (10.2) 

Although one has B • J = 0, nothing can be deduced however for the scalar 
product B • V. 

By vector multiplication of (|9.9|) by V, we get: 

Vx— = JvxE - c 2 B + (V-B)V) (10.3) 



Dt 

that generalizes ()6.19j) . Finally, by scalar multiplication of l|9.4[) by E and 
1)9. 5 j) by B, one obtains: 

1 d 

- — (|E| 2 + c 2 |B| 2 ) = -c 2 div(ExB) - c 2 (V-B)divB (10.4) 

which is the counterpart of (|3.5I) . 

Referring to figure 7, let J^^ be the normalized Poynting vector at time 
t and j( t+5t ) the one at time t + 5t. Let V be the vector at time t, obtained 
by backward parallel transport along the stream-lines of the vector 
Then, we have: 

DV V-cjW 

— - = lim 10.5 

Dt 6t-*o dt K ' 

We recall that the same was done in section 7 in order to define the vector G 
(see (|7.3[) ). Therefore, for small time variations St, we are allowed to write: 

V Rs cJ + G5t = cJ + n(E + cJ x B)5t (10.6) 

with J = (E x B)/|E x B|. If E • B = 0, the relation (fTTJTE|) can be rewritten 
as: 

V ss cJ + fj, (|E| - c|B|)(Jt (10.7) 
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Figure 7: Typical behavior of different fields, when a wave-front is forced to turn. 

after noting that: (E x B) x B = (E • B)B - |B| 2 E = — |B| 2 E. This shows 
that it is sufficient to have |E| 7^ c|B|, in order to activate the transversal 
field G. 

We can compare the evolution of an electromagnetic phenomenon to 
that of an inviscid fluid, whose mass density, up to dimensional constants, is 
given by p = divE. Note, however, that a real "mass" does not exist. Note 
also that p can also attain negative values. We do not define the density 
p = divB for reasons that will be detailed at the end of section 14. The 
following continuity equation holds (see also (|3.4j0 : 

^ = -div(pV) (10.8) 

obtainable by taking the divergence of (|9.4jl . The equation l|10.8j) can be 
also written as: 

^ = -pdivV (10.9) 

Dt H y ' 

For a plane wave (or soliton) having p ^ 0, we obtain divV = as well as 
G = 0. Then, (|10,9[) tells us that the fluid shifts, without modifications, 
along the direction determined by V. The fluid travels at constant speed 
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c, showing rarefactions and compressions,. More properly, it evolves like 
an incompressible fluid, but with density not equally distributed in space. 
Regarding a spherical wave, one has divV > and G = 0. As expected, 
this implies that the mass density diminishes (in absolute value), while time 
passes, since it spreads on spheres of growing area. In both examples (the 
plane and the spherical) we have curlV = 0. In other words, the fluid is 
irrotational. 

Let us suppose that V is a gradient, i.e.: V = V\&, where ^ is a scalar 
potential not depending on time. Then, the corresponding fluid is irrota- 
tional. Thanks to (|9.7jl . we trivially have: 



which is the eikonal equation. Then, we observe that (|l(J.lUj) and the relation 
4* = 0, imply: 



This confirms a remarkable result: the eikonal equation (hence, the evolution 
of the wave-fronts based on the Huygens principle) is perfectly compatible 
with the condition G = 0. This analytic property, obtained without approx- 
imations, goes beyond the famous limits of geometrical optics. So, here, with 
a very elementary proof, we obtained another important result. 

In the equation (|9.9j) . the term VxB = — BxV = T(V) can be 
viewed as a suitable stress tensor T applied to the vector normal to the 
front of propagation (see for instance JT|, p. 10). As we already know, forced 
variations of the electric field produce changes in the motion of the fronts. If 
these are combined with forced variations of the magnetic field, a torsion is 
also introduced, which modifies the polarization of the wave. We guess, that, 
when the external perturbations end, the electromagnetic fields return to 
their natural equilibrium in which |E| = c|B| and (|6.18|) is satisfied, so that 
the fluid takes again an irrotational motion. From the examples discussed in 
sections 7 and 8, this behavior corresponds to what is commonly observed 
in nature, and certainly comes from the minimization of some Lagrangian. 
At the moment, we do not however have theoretical explanation for this 
conjecture. 

Furthermore, we note that stationary electric fields, for example with 
B = 0, are not longer solutions. We can check this by examining relations 
(|9.9|) and (|1U.2|) . They force the velocity field V to turn itself around (V 



W| = c 
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deviates in the direction of E, but E remains orthogonal to it). More gener- 
ally, equation Q9.6JI requires the solutions to be in continuous evolution. We 
contend that the new system of equations does not admit stationary solu- 
tions having finite energy. We made the same consideration in section 5, in 
the particular case of solitons. Nevertheless, there could be nonstationary 
solutions localized in space. We can imagine for instance the case of two (or 
more) solitons, in such a way that they are constrained, by influencing their 
electromagnetic fields each other, to revolve around a common center. We 
still do not have all the elements to study these phenomena, which, as we 
will see in the following pages, need the environment of general relativity to 
be stated properly. Some remarks about the case of two rotating solitons 
will be given in section 15. 

We are unable at the moment to obtain the equations (|9.4JI - (j9.5JI - (j9.6j) 
from the minimization of a suitable action function as we did in section 6 
(concerning (j9.fi j) alone, something in this direction will be obtained in the 
next section). One may consider the usual Lagrangian L = 2(c 2 |B| 2 — |E| 2 ) 
and the generalization of the relation (|dl2[) . i.e.: 



Then, we could minimize the same action function given in ()6.6|) using the 
constraint ()1(J.12|) . Nevertheless, we would not obtain the desired result, 
since (jlU.12|) is too restrictive. In this way, we only get a set of equations 
describing free waves. As a matter of fact, we can prove that, if V has 
the same direction as A, then one automatically has G = 0. Assuming 
divB = 0, this check can be done as follows. Considering (|6,1|) and Q1U.12JI . 
we have: 






M V<f + i(A-V)A 



Dt 





DV 



) 



DV 



/z(E + VxB) = (10.13) 



Dt 



- V 



E + 



—Ax curl A 



Dt 
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where we used that <fr 2 = |A| 2 . The last equality is true thanks to (|9,9|) . 
Then, noting that (c/<& + p,/c)A = (1 + ^<£/c 2 )V, the relation ()10.13|) leads 
to: 



D 
Dt 



V + 1 + £=- — = (10.14) 



c 2 Dt V c 2 J Dt 

By scalar multiplication by V, due to (|1U.1|) . we recover: 

$£"» + (. + £)£-v-,£-. (-) 

Thus, $ turns out to be constant along the stream-lines. For this reason, 
from (|1L). 13|) . also A is constant along the stream-lines. Therefore, one has 
-j^V = G = 0. We also conclude that, when the rays bend (G / 0), then 
the vector A cannot be aligned in the direction of motion. 

We mentioned in the previous sections that the mechanical effects are 
implicitly included in the term c 2 divP, where it is necessary to distinguish 
between the contribution due to the variation of the Poynting vector along 
the actual direction of propagation of the front, and the transversal contri- 
bution (which is zero when G = 0). Differentiating with respect to time the 
expression J = P/|P|, we get: 



G = 1ST -5T " ^T^ P ( 10 - 16 ) 





In particular, by scalar multiplication of G by P, (|10.16j) shows the orthog- 
onality relation G • P = 0. Furthermore, from (|H.5|) . the energy can be 
described as a work by integrating — 2c 2 divP with respect to time. This 
yields: 

- 2c 2 J divP dt = -2 J divP V • ds (10.17) 

where we set ds = \~dt, which implies V • ds = |V| 2 cit = <?dt. The last 
integration is made along the curve T representing the path of the light ray. 

We end this section by illustrating another interesting relation. Let us 
define as usual p = divE. Afterwords, let us assume that p ^ and define 
lD = F/p, where F = curlV + /iB. Then, along the stream-lines we have: 

^ = (B-V)V (10.18) 

Note that u> is dimensionally equivalent to a time multiplied by a charge 
and divided by a mass. The equation (|l().18j) recalls the analogous one for 



45 



isentropic flows, which is introduced in fluid dynamics by defining u> as the 
curl of the velocity field divided by the mass density (see [5], p. 24). Using 
()6.1|) . the field Q also takes the following form: 



curl(V + /uA/c) 



ia_ 

cdt 



(10.19) 



divA - A$ 



The equation (|1().18|) can be proven as follows: 



§-(,-v,v 



1 DF 

~p~Dt 



J_£p f 

P 2 Dt 



(F • V)V 



curl 



f ~Qf ) + + (V • V)F + F divV - (F • V)V 



dt 



curl( - (V • V)V + p(E + V x B)J + p(- curlE - V divB) 

+ (V • V)F + F divV - (F • V)V 
= curl(V x F) - V divF + (V • V)F + F divV - (F • V)V 

- curl[(V • V)V] + curl(V x curlV) 

1 



curl((V • V)V + (V x curlV) 



curl(V|V| 



(10.20) 



where, in the order, we used (|10.9|) . (|9.6() . (|9.5|) . some well-known calculus 
properties and the fact that V|V| 2 = 0. 

In the case of plane solitary waves, we have curlV = 0, hence <D = pB/p 
(when p ^ 0). Therefore, u remains orthogonal to V, so that the relation 
(|10.18j) becomes -j^-cj = 0. Then, the quantity u shifts, remaining constant 
along the stream-lines determined by the velocity field V (which are straight- 
lines in this case). 



11 Towards general relativity 

Our first step, in this section, is to recover the equation (|9.9|) through the 
minimization of a suitable Lagrangian. To this end we work in space-time 
using 4-vectors. Let us start by defining (xq, x\, X2, 23) = (ct, —x, —y, —z) 
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and (e , ei, e 2 , e 3 ) = (1,-1,-1,-1). Then, for the vector (V , V\, V 2 , V3) = 
(Vo, V), one has: 

3 

5>lf = V 2 - |V| 2 (11.1) 

i=0 

As in section 6, we assume that divB = and introduce the potentials <3? 
and A by q6.1|) . Let also be (A , A\, A 2 , A 3 ) = ($, A). Up to multiplicative 
constants, we can define a Lagrangian in the following way (see also [TU] . 
p.50): 

L = c v^o 2 - |V| 2 + n ($ - i A • (11.2) 

The quantities Vi, i = 0,1,2,3, are the independent variables, while the 
potentials depend on Xi, i = 0,1,2,3. By setting Vq = c, the term in 
parenthesis of (jll.2f) can be written as: c" 1 Yh=q ^i-^-iVi- For the moment, 
we do not impose the condition (|9.7|) . implying that the sum in is 
zero. 

Suppose that we are moving along a stream- line (or curved ray), between 
two instants of time t± and t 2 , then the action function takes the form: 

S = - f t2 Ldt (11.3) 



d ( dL\ dL dL d ( dL\ dL . A . 

-(-)=^ = c ^ dt{dvj = -d^ * = 1 ' 2 ' 3 (1L4) 

where we observed that -V) = (^, In particular, for 



d ( dL \ d ( —cVi fi 

— — A; 



Its minimization brings to the Euler-Lagrange equation (see jH], p. 577): 

d ( dL\ _ dL _ dL d / dL\ dL 

M\d%)~~dt~ C dx^ Jt\dVi) ~ ~~dxi 

! ( 

i = 1,2, 3, we have: 
d_ 

dt \dViJ 

\ y ' u i 

DVj _ vdA^ _ /Xy^SAj dxk_ (115) 
Dt c dt c ^ dx k dt ' ' 

dL d$ (i dA 

and — = (j,— - - V ■ — (11.6) 

oxi dxi c dxi 

where, in (|11.5|) . the substantial derivative -SiVi gives the variation, along 
the stream-lines, of the coordinates of the velocity field, parametrized with 

1 f* . A/2 _ IVI2, 




respect to the arc-length: s = \c\ J t JVq — |V| 2 (i£ 
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If we now define = for k = 1, 2, 3, thanks to (|6.3|) and (|6.5j) . we 

get: 

DVi ( 1 9 A 



c <9t 

- H V • fvA + ^) = - £ (11-7) 

When Vo = c, the last term in (|11.7|) is equal to the i-th component of the 
vector /i(E + V x B). This implies the equation ()9.9|) , 

Concerning k = 0, if we fix Vq to be constantly equal to c, one obtains 
-jj[Vo = 0. Therefore, we have: 

d ( 8L\ DV d$ <9$ „ , 1iB , 

+ A 4 ~rr = "rrr + /xV-V$ (11.8) 



\dV J Dt r dt r dt 

dL d$ u „ <9A , . 

Due to 1)11.4(1 . by equating these two last expressions, we recover: 

= -ufv* + -^Vv = = ^E-V (11.10) 

\ c at J 

which corresponds to (|10.2j) . Considering ()11.10|) . by scalar multiplication of 
(|9.9|) by V, we deduce that the field minimizing the action satisfies V-jpgV = 
0. Hence, the norm | V| is constant. If such a constant is c, we finally obtain 
the relation (|9.7jl . that says that the solutions evolve on the light-cone. 

At this point, it should be noted that, by choosing |V| 2 = c 2 , the first 
part of the Lagrangian in (|11.2|) vanishes. This does not mean that it van- 
ishes identically, but only in correspondence to the minimum. Instead, the 
second part of the Lagrangian is zero when A- V = c$, which is very similar 
to the condition (|6.17|) . obtained from the constraint (|6.12|) (see also (|1U.12[0 . 
This coincidence is quite significant. Perhaps, in the future, it will suggest 
a way to build a Lagrangian for the entire set of equations (|9.4|) - (|9.5j) - (|9.6[) . 

By multiplying the equation (|11,7|) by Vi, i = 1,2,3, and the equation 
(HTTPl) by V , we get: 

F ik V k Vi = (11.11) 

where the sum is for i and k going from to 3. This also trivially follows from 
the anti-symmetry of the tensor F . The equation (|dll|) is also written as: 

(Qpik \ / QJpOk \ 

~dx~) V ° ~ \~dxTj 6lVi = ° to * = (11-12) 
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Otherwise, the equations ()2.3|) and (|2.4|) . can be recovered from the expres- 
sion (see for instance UJ, p. 150): 

„ uFjk oFu-i oF™ , . 

F ikj = + -k^ + -tr 1 = 11.13 

oxi ox k 

where there is no sum on repeated indices. The rank-three tensor F^j is 
anti-symmetric and called the cyclic derivative of Fn~. On the other hand, 
the equation (J9.5|) follows from the expression: 

'dFi k , dF kj dFjA _ fdF 23 dF 31 dF 12 



vo + + ± + + 

\ axj axj OTfc / V ctei ax2 0x3 / 

where the indices m,j,i,k (taken in this order) are all different. The sign 
± depends on the permutation (even or odd) of the indices (the sign is plus 
if m = 0, j = 1, i = 2, k = 3). In (|11.14j) the term in parenthesis on 
the right-hand side is equal to cdivB. In a more contracted form, the last 
equation reads as follows: 

Vb F jik = ± e m y m F123 (11.15) 

In the results obtained above, we basically considered V as the velocity 
field of an infinitesimal particle moving at the speed of light. On the other 
hand, in a wave there are infinite contiguous trajectories. As a matter 
of fact, the evolution of a wave is a global phenomenon, that should be 
taken as a whole, and not studied independently along each stream-line. 
For such a more in depth analysis, we need to work in the context of general 
relativity. We are going to show that the passage of a wave modifies the 
space-time structure. For a free wave this does not affect the evolution of 
the wave itself (see section 13), but for constrained waves the change of the 
geometry influences their entire behavior. The analysis will allow us to find 
the coupling between the fields describing the wave and space-time geometry, 
hence the link between electromagnetic and gravitational phenomena. 

We first need to introduce some classical definitions. Mainly, we adopt 
the notation used in [Jj . The space-time geometry is locally determined by a 
symmetric bilinear form (the metric tensor), whose coefficients are denoted 
by gij. Then, the Christoffel symbols are defined in the following way: 

F i _ g im (dgmk dg mj dg kj \ 

" —1^7 + ^ ~ dx~ (1L16) 
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where we sum over the index m. The coefficients g %3 are in such a way that: 

g im g mj = Sij (11.17) 

The coefficients are adimensional, while the Christoffel symbols are the 
inverse of a distance. If the space is "flat" (Euclidean or Minkowski space), 
all the Christoffel symbols vanish. In this case, one has g tk = gik = eiSik- 
As usual, we denote by g the determinant (which is negative) of the tensor 
gik- A lemma due by Ricci (see j7], p. 129) claims that the 4-divergence of 
the metric tensor is zero. In detail, one has: 



v ,„ I d{^-gg^) + n ( j i ]S) 

-g dx k 



jk _ 1 u w~yy ) , r i Am 



where is the covariant differentiation operator. The same is true for the 
coefficients g^. Moreover, the coefficients g lk are said to be harmonic when: 

v ' ■' = (11.19) 



V-g dxk 

Next, we define V 1 = g tm V m . The values V m are the entries of a velocity 
vector expressed in the coordinates system xo,x%, X2,x%. Let us set V° = c. 
Then, the condition (|9.7|) is generalized in the following way: 

V% = g ik V k Vi = gimV'V™ = (11.20) 

In this more general framework, the equations (|11.7|) and (|ll,10j) are 
rewritten as: 

^ + Y) k V 3 V k = -^F™V m fori = 0,1, 2, 3 (11.21) 
JJt J c 

For i = 0, 1, 2, 3, we also define (see [7j, p. 217): 

G i = ^^ + T i jk V j V k = V m V rn V l (11.22) 

From (|11.21|) we easily recover the orthogonality relations: 

GiV* = G% = (11.23) 

where Gi = gi m G rn . Finally, let us define G = (G\, G2,G^), which is dimen- 
sionally equivalent to an acceleration. 
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In general relativity, the gravitational field is somehow identified with 
the tensor gij. Of course, the vector G may vanish without having that 
the space is flat. Although G does not fully characterize the properties of 
space-time, for us it will be the "real" gravitational field, that is the one 
we can measure in our everyday life. In the following, G will be called the 
vector gravitational field. 

The equation (|11.21j) enables us to understand how the trajectory of a 
"thin" solitary wave can be distorted when immerged in a given gravitational 
field (having G / 0). Being the soliton a free wave, the right-hand side of 
(|11.21j) vanishes (see (|6.18|) V Thus, its path follows a suitable geodesic in 
space-time, the shape of which is determined by the external gravitational 
field. This should correspond to some transversal bending in the direction 
locally individuated by the vector G. With this reasoning, we have to ne- 
glect a couple of facts, both due to the change of direction: the modification 
of the electromagnetic fields and the "gravitational reaction" of the soliton 
(a curving wave produces new gravitational field). As we argued in sec- 
tion 10, these should be minor effects, since the wave, for some principle 
of least action, tries to compensate the electromagnetic fields in order to 
enforce (|6.18|) . From the point of view of the soliton, the path followed is 
straight, even if it actually travels on a curved geodesic. To get more reliable 
quantitative results, we must solve a quite complicated system of equations. 
Comparing the computed results with the experimental evidence, one could 
probably evaluate the constant /x. This is however an exercise that we would 
prefer to avoid here. Also if some theoretical passages may be formally sim- 
ilar, the important clue is that there is no need to suppose that a soliton 
has an infinitesimal mass to justify that is attracted by a gravitational field. 

After recalling that F tk is an anti-symmetric tensor, in general coordi- 
nates, the equation (|11.12j) becomes: 

■te^.^f*!^, i = 0,1,2,3 (11.24) 
V=5 V ox k dx k J 

or, in more contracted form: 

{V k F ik )V° - {VkF^V 1 = » = 0,1,2,3 (11.25) 

The equation (jll.LSj) remains unchanged. However, it can be also written 
in the following way (see [7j, p. 133): 

F lkj = VjF ik + V t F kj + V k Fji = (11.26) 
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Besides, equation (|11.15|) becomes: 

V° F jik = ±V m F 123 (11.27) 

By taking the 4-divergence of the contravariant vector in (jll.24j) and con- 
sidering once again that the tensor F lk is anti-symmetric, we arrive at the 
continuity equation: 

i *(,/^»V)_ w , th pE i gu^j 



We got an analogous result in section 3, by taking the standard divergence of 
the vector equation (|3.1|) . The time derivative came from the term div^J^-E^ 
and the term div(curlB) was zero. All the pieces here combine in a com- 
pletely different manner. Nevertheless, the final result is extraordinarily 
similar. 



12 The energy tensor 

Let us first work in Cartesian coordinates. We will define the symmetric 
electromagnetic stress tensor in the classical way (see [7j, p. 96), i.e.: 

Utk = -(y^ejFijFkj - i(c 2 |B| 2 -|E| 2 ) ej 5 ife j (12.1) 

We have U 00 = |(|E| 2 + c 2 |B| 2 ) and £f=o BiUu = 0. 

Its contravariant version is given by U lk = e^e^U^ and of course we have 
J2i=o e iU n = 0. The explicit expression of the contravariant tensor is the 
following: 

/i(|E| 2 +c 2 |B| 2 ) cB 2 E 3 -cE 2 B 3 cE^-cB^ cB 1 E 2 -cE 1 B 2 \ 
cB 2 E 3 ~cE 2 B 3 lf ( + 2 f B f-| C EP) ~E l E 2 -c 2 B l B 2 -E X E 3 - c 2 B x B 5 
cE 1 B 3 -cB 1 E 3 -E 1 E 2 -c 2 B 1 B 2 lf ( + 2 f B f ~E 2 E 3 -c 2 B 2 B 3 
\cB 1 E 2 -cE l B 2 -E x E 3 -c 2 B x B 3 -E 2 E 3 -c 2 B 2 B 3 ~_f£ f B f j 

If ()9.4j) and ()9.5j) are satisfied, then an important property of this last 
tensor is that, in the case of free waves (hence in absence of mechanical 
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terms), its 4-divergence vanishes. Indeed, we have for i = 0, 1, 2, 3: 



QJjik 

provided (|6.18|) (or (|6.19j) ) is satisfied. This implies that E • B = and 
|E| = c|B|. These hypotheses also imply that |ExB| = |E||B| and V-B = 0. 
Let us prove (|12.2fl starting from i = 0. Thanks to (j3.5fl . one has: 

BTJ ok 1 r) 

fcT = S» ( l E|+c l B l ) + cdiv(ExB) = (12 ' 3) 

As far as the other values of i are concerned, let us begin to define: 

N = (Ni,N 2 ,N 3 ) = — c 2 curlB + (divE)V 

ot 

M = (M 1 ,M 2 ,M 3 ) = — + curlE + (divB)V (12.4) 

Thus, if the equations (|9.4|) and (|9.5|) are true, then we get N = and 
M = 0. We are ready to check (|12.2|) for i = 1 (the other cases are treated 
in a very similar way). We have: 

dU lk d , d 



Qxk dt {B 2 Ez-E 2 B s) --(-Et + ^ + c*Bi) 

+ \ ^( C2|B|2 " |E|2 ) " T { ~ ElE2 ~ ° 2BlB2) " §- Z { ~ ElE3 ~ ° 2BlB3) 

= (M 2 E 3 - M 3 E 2 + N 3 B 2 - N 2 B 3 ) + (Ex + V 2 B 3 - V 3 B 2 ) divE 

+ (c 2 Bx + V 3 E 2 - V 2 E 3 ) divB = (12.5) 

The last three terms in (|12.5p are actually zero for the following reasons. In 
the first one we recognize the second and the third components of N and M, 
which vanish, if we assume that the equations ()9.4|) and ()9.5j) are satisfied. 
The second one contains the first component of the vector E + V x B, which 
vanishes in the case of a free wave. Concerning the last term, the part in 
parenthesis is the first component of the vector c 2 B — V x E, which also 
vanishes (see (|H.19jl ). 

The property (|12.2|) is reported in many texts (see for instance T , p. 97). 
But, it is extremely important to observe that, in the case of Maxwell equa- 
tions, the last two terms are zero because divE = and divB = 0. This is 
the reason why we decided to double check equation (|12.2j) , which turns out 
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to be fulfilled even when the divergence of the fields E and B is not zero 
(the assumption we supported throughout this paper). Therefore (|12.2|) 
holds under weaker hypotheses. 

As expected, a converse statement also holds: assuming that (|12.2|) is 
true, then we can recover the equations (|9,4|) and (|9.5|) , This amounts to 
differentiate the equation of energy conservation, in order to obtain the 
corresponding Euler equations. Arguing as we did to get (|12.5|) . we arrive 
at: 

/ dU lk dU 2k dU 3k \ 
\ dx k ' dx k ' dx k J 

+ (E + V x B) divE + (c 2 B - V x E) divB (12.6) 

Assuming, as previously done, that we are dealing with a free wave, after 
eliminating in (|12.6j) the vanishing terms, we are left with (M x E — N x B). 
Since, by hypothesis, the equation (|12.2|) is true, if the vector N is zero, 
then M must be also zero (likewise, if M is zero, then N is zero). Therefore, 
(|9.4I) is satisfied if and only if Q9.5|) is satisfied. This is the same situation 
encountered in the classical Maxwell equations, where J^B + curlE and divB 
both vanish if and only if J^E — c 2 curlB and divE are both zero. In the 
standard approach, the first pair of equations are satisfied by choosing the 
pontentials A and <3? as in ()6.1|) . The second pair is obtained by means of 
variational type arguments. 

Of course, we can find "intermediate" situations, by suitably redefining 
the two potentials. Let us take for example: 

1 OA 

B = curlA £ = V$ 

c dt 

= AE + (l -A)cB = AcB - (1 - A)E 

^ + (1-A) 2 ^X' + (1 - A) 2 

where A is a real parameter. From the relations (|12.7|) we can explicitly 
compute the fields E and B in terms of £ and B. These also imply: 

dB 

div8 = — = -ccurLS (12.8) 

dt y ' 

Then, it is a matter of minimizing the usual Lagrangian. At this point, 
introducing the constraint A • V = <£, one gets the equation: 

3£ 

— = c curLB - V div£ (12.9) 
dt 
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that, for A = 1, is equivalent to equation (|9,4|) . The equations in Q12.8JI are 
equivalent to require (A — 1)N + AM = 0, while the one in (|12.9|) brings to 
AN + (1 - A)M = 0. 

It is to be noted that V has the same direction of E x B, which is also 
like that of £ X B. So, from the energy tensor it is not possible to figure 
out what the parameter A is, as well as the polarization of the wave. This 
information has to be provided with the initial conditions. For instance, the 
wave in (|7.4|) . circolarly polarized, produces the same tensor U lk of a linearly 
polarized wave, moving in the same direction with twice the intensity. As 
a further consequence, we finally observe that U tk does not change if E 
takes the place of — cB and cB takes the place of E. Such a permutation 
corresponds to the choice A = 0. 

We can now argue in a general framework. For a given metric tensor go- , 
the electromagnetic stress tensors must be modified in the following way 
(see 0, p.151): 

U ik = - (g mj F im F k i - \g ik F mj F^) (12.10) 
where F^ is given in (|6,4j) . while F lk comes from the relation: 

F ik = g im g kl Fmi (12 n) 

Assuming to be as in the case of a free wave, the equation (|12.2j) has to be 
replaced by the following one: 

V k U ik = - W / ' + T l m AJ mj = (12.12) 

V^g ox k J 

The proof of (|12.12|) is given for instance in 9[, p. 606, for the classical 
Maxwell equations. This is also true in the case of our new set of equa- 
tions (hence under weaker hypotheses). At the end of this section, we will 
evaluate the 4-divergence of the tensor U lk in general coordinates. Such 
generalizations are unavoidable since, even the simple case of a plane wave 
provokes a modification of the space-time geometry, requiring to work with 
tensors of the form (|12.1()|) . These aspects will be better studied in the next 
section. 

When the electromagnetic phenomenon is not a free wave, we cannot 
expect that (|12.2|) and its generalization (|12.12j) are verified. This means 
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that the system constituted by the sole electromagnetic part is not energy 
preserving. We know that, in this case, the energy balance has to take care 
of the mechanical effects. Thus, we study how to introduce them. We start 
by assuming that divB = 0, leaving the discussion of the more general case 
to section 14. Then, let us define a mass tensor as follows: 



The contravariant version is given by M 
written as: 



M ik = V l V k divE (12.13) 
— ejefcMjfc, which is explicitly 



M 



ik 



V 2 


-V Vi 


-v v 2 


-V V 3 \ 






v x v 2 


V X V 3 


V V 2 


V1V2 


v? 


v 2 v 3 


v v 3 


VtV 3 


v 2 v 3 


vi J 



(12.14) 



where Vq = c and p E = divE is a kind of mass density (dimensionally this 
is not correct, but this aspect will be altered later). We recall that p E can 
also be negative. Let us check what happens to gf-M l ' fc . For i = we have: 



8M 0k 
dx k 



c c 



dpE d{p E Vi) d(p E V 2 ) d(p E V 3 ) 



8xq dx\ 
dp E 



dx<2 



dx 3 



at 



+ div(p E V)^) = 



(12.15) 



This is true because of the continuity equation (|1U.8JI with p = p E . For the 
other indices i = 1,2,3, we have: 



g M ik 
dx k 



d{ PE Vi) d(p E ViV i ) d(p E V 2 Vi) 



d( PE V 3 Vi) 



dx 



dx\ 



( d( PE Vi) 

\ at 



PE 



+ div(p E ^V) 



dt 



+ (y-v)Vi 



dx 2 

dp E 

dt 

DVi 
Dt 



dx 3 

+ div(p E V) 
divE 



(12.16) 



where we again used the continuity equation. We conclude for instance that, 
if the light rays are straight-lines (that is: G = 7^7 V = 0), then one gets 
9-M ik = 0, for i = 0,1,2,3. 
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In non Euclidean geometry, it is necessary to generalize the mass tensors 
in the following way: 



1 F ok ) 



M ik = pvViV k M ik = p- E V i V k with p E = ' " vv ' (12.17) 

V 9 ox k 

With the help of the continuity equation 1)11.28(1 and the definition 1)11.22(1 . 
it is easy to get, for i = 0, 1, 2, 3: 



1 d(^g M ik ) 



dx k 



PvV k j— + pA^V* = peG* = - ^ PE F ik V k (12.18) 
ox k J c 



where the last equality is a consequence of 1)11. 21|) . Hence, in a flat space 
{G l = 0), the 4-divergence of the mass tensor vanishes. Moreover, we observe 
that the mass tensor does not contain the pressure term (on the other hand, 
an equation of state is not defined). 

In order to combine electromagnetic and mechanical effects, we sum up 
the corresponding tensors, by defining: 

T ik = ^(pU ik + M ik ) (12.19) 

where the constant p is dimensionally equivalent to a charge divided by a 
mass. It follows that Ti k has the same dimension of a curvature, that is the 
inverse of the square of a distance. Now, in a flat space-time, even if we are 
not dealing with a free wave, we may write: 

— - = (12.20) 

dx k 

As a matter of fact, due to 1)9.9(1 . if the term p(E + V x B)divE of the 
electromagnetic part does not vanish (see 1)12.6(1 ). it is anyway compensated 
by the corresponding term — divE-^V of the mechanical part (see (jl2.1fij) ). 

In the general case, the relation ((12.2UI) is substituted by: 

V k T ik = VV y '- + r^ T"^' = (12.21) 



Before ending this section, we would like to verify that ()12.21)1 actually 
corresponds to the Euler equations. As a matter of fact, ()12.21|) is satisfied 
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when pi.25|) . (|11.27j) and Q11.21JI are true. We recall that these three last 
equations are the generalizations of (|9.4|) . (|9.5|) . (|9,6|) . respectively. For the 
moment, we will only treat the case in which divB = 0, leaving the general 
discussion to section 14. Let us start by computing the 4-divergence of the 
tensor U lk . First of all, we have: 

V k U lk = -V k {g mj F im F k i) + \V k {cf k F mj F m i) 

= 9 m V k (F mi Fi k ) + \g lk V k {F mj F m i) (12.22) 

where we notice that g mj F im = g im F mj (thanks to (112. Ill) ), that F k i = 
-Fi k and that, due to (111. 1811 . it is allowed to exchange the metric tensor 
with the covariant derivative (see also [7j, pl30). Going ahead, one gets: 

V k U ik = g m (V k F^ k )F mj + g m (V k F mj )F^ k + y k {V k F m] )F m ^ 

= c-\V k F» k )F m3 g m Vi + \g im {V k F m3 )F^ k 

+ \g im (V k F m] +V m F 3k )Fi k (12.23) 

where we used (lll.25|) (with V° = c). The other passages have been obtained 
by a suitable renaming of the indices. Recalling the definition of Fi k j given 
in H11.26|) . we have: 

V k U ik = ^F m3 g m Vi + y m (V k F mj )F^ k 

+ y m F m]k P k - y m (v,F km )^ k 

= ^F im V m + \g im F m3k Fi k (12.24) 

In the last passage two terms have been deleted, since, after renaming the 
indices, they resulted in being equal and with opposite signs. The last term 
in (|12.24[) is zero because of (|11.2fi|) (remember that we are studying the 
case divB = 0, thus -F123 = 0). Of course, the final result is zero when, for 
instance, pe = 0, as in the classical Maxwell case. But it is also zero when 
F im V m = 0, which corresponds to the case of a free electromagnetic wave. 
On the contrary, we need to consider the contribution of the mass tensor. If 
F m j k = 0, taking into account the relations (|12.19|) and (|12.18|) . one finally 
obtains: 

V k T ik = £ (pV k U ik + V k M ik ) 

= 4 (-PEF im V m ~ ^ PE F* k V k ) = 
c 4 V c c / 
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13 Unified fields equations 



In the previous section, we build the symmetric tensor that includes both 
the energy contribution of an electromagnetic wave and that of mechanical 
type, taking into account possible deviations from the natural propagation 
path of the wave. The properties of insure the preservation of energy and 
momentum. Hence, we can put on the right-hand side of the Einstein 
equation: 

Rik ~ \gikR = X Tik (13.1) 
in which we recognize the Ricci tensor: 



Rik 



dx k 



L ik L jr 



-pj pm 

im kj 



(13.2) 



the scalar curvature: 



R 



g lk Rik 



113-3) 



and an adimensional constant x- We recall that the Christoffel symbols 
are defined in (|11.16|) . Let us note that, by (|12.10f) and (|12.17j) . the metric 
tensor, which is now our unknown, also appears on the right-hand side of 

jmj). 

We soon examine the response of equation 1|) to the passage of the 
most elementary plane wave. We take for instance the expression given in 
(|2.5jl . where we have E\ = cBi = csinu;(i — z/c), divE = and divB = 0. 
We will verify that, even in this simple case, the space-time geometry, that 
comes from the solution of ()13.1|) . is not Euclidean. In fact, we look for a 
metric tensor g^ of the following type: 



gik 













-p 2 











-1 


Vo 









J,k 











°\ 






-l/p 2 











-1 





Vo 








-1/ 



;i3.4) 



where p is a function, to be determined, of the variable £ = t — z/c. Some- 
how, we are expressing a preference for the direction of the x-axis, which is 
orientated with the electric field. The determinant g of is equal to —p 2 . 
The corresponding Christoffel symbols are: 



r o _ PP 



r 3 



-pp' 



ii 



01 



10 



13 



r i 

isi " cp 



;i3.5) 
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where the prime denotes the derivative with respect to £. All the other 
symbols vanish. The non-zero coefficients of the Ricci tensor are instead: 



R 



oo 



R 



03 



R 



30 



^33 — n~ 

c L p 



;i3.e) 



The scalar curvature R is zero. 

Being zero the divergence of E, the mass tensors e M tk vanish. 
Actually, one should check that pe = 0. This is also true, as the comments 
at the end of this section illustrate. The tensors and U lk have to be 
computed through (|12.10|) . First of all, one has: 



F 



4k 











-u 










u 
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Fik = 


c 

















V0 


-u 


0/ 
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u/p 2 









u/p 2 
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—u/p 2 






113-7) 



:i3. 



/ 



where u = B 2 = E x /c. Note that (V , V) = (c, 0, 0, c) and (V°, V 1 , V 2 , V 3 ) = 
(c, 0, 0, — c). Hence, one gets F tk Vk = 0, from which we deduce that the wave 
is free, as is already known. Afterwards, we have: 



T. -fL 

- L ik — o 



/ (u/P? 




\(u/p) 2 















\-(u/pf 













{u/p) 2 \ 



(u/p) 2 J 

-(u/p) 2 \ 



(u/p) 2 



(13.9) 



(13.10) 



Thus, Q13.1JI and ()13.10|) bring us to the equation: 



p"p 



2 2 



(13.11) 



For u = smuj(t — z/c), we finally obtainp = (u^/x/u) sina^i — z/c), which is 
the solution we were looking for. There are surely other geometries compati- 
ble with the same plane wave. Note that the one presented here satisfies the 
relation (|11.19|h We also observe that there are points in which the metric 
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becomes singular, that is, the determinant g is zero. Another equivalent 
possibility is to exchange gu and #22 m (|13.4|) . and make p oscillate with 
the magnetic field. Comments about this option will be given in the next 
section. 

The solution just obtained can be assimilated to a transversal (perfectly 
plane) gravitational wave, travelling in phase with the electromagnetic one. 
It must also be noted that, even if the space is officially non Euclidean, the 
geodesies involved in the motion of the wave are straight-lines. The field 
G, defined by (jll.22j) . is identically zero. This is in agreement with our 
viewpoint: the geometry may be deformed, but there is no creation of a real 
gravitational vector field. 

Pure gravitational solutions resembling plane waves, were formerly de- 
tected in j3] . We have been able to get the above explicit (and very simple) 
solution because we were resolute enough to assume the dependence from 
the metric tensor of the right-hand side of the Einstein equation. As far as 
we could deduce from the current literature, in contrast to our general ap- 
proach to the problem, it is customary to construct the electromagnetic en- 
ergy tensor in vacuum (thus, in Minkowski space-time), also because such an 
assumption is supposed (erroneously) to simplify the computations. Then, 
one comes to a set of solutions, but, as we proved, this is not the correct set- 
ting. Note also that, commonly, gravitational waves are searched among the 
solutions of the linearized homogeneous Einstein equation, obtained after 
perturbation of the flat space-time. 

We can recover the laws of motion by evaluating the 4-divergence of the 
tensor T lk in (|13. 1(J|) . The geometry is non Euclidean, therefore, the relation 
(|12.2()j) has to be replaced by (|12.21[) . where \f—g = \p\. For i = 0, one has: 



1 d(^=g~ T ok ) 



-(J 



dxi 



+ 



pU rpm] 

mj- L 




c 2 p 



The situation is exactly the same for i 
when for instance: 







c 2 p 






f 1 dp 


p2 1 


K c dt 



at 



+ 



+ 



d(u 2 /p) 
dz 



dp 

dz 



(13.12) 

3. The last term in (113.1211 is zero, 



1 du du q 
c dt dz 



1 dp dp 
-— + — = 

c dt dz 



(13.13) 



Our plane electromagnetic-gravitational wave is actually the solution to both 
the above equations, once the proper initial conditions have been assigned. 
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Let us now discuss the example of a circularly-polarized plane wave: 



E = (ccos tu(t — z/c), csinuj(t — z/c), 0) 
B = (— sinuj(t — z/c), cos u(t — z/c), 0) 



(13.14) 



The classical divergence of the electric field vanishes, as well as that of the 
magnetic field. Then, let us take the following metric tensor: 



9ik 



M 2 X 

LU 2 



/I 




Vo 





[cos u(t — z/c)] 2 








[s'muj(t - z/c)] 2 




° \ 




-1/ 



(13.15) 



in such a way that the coordinates x and y are syncronized with the electric 
field (the reasons for this choice will be explained at the end of section 14). 
In this case, the nonvanishing coefficients of the Ricci tensor are: i?oo = 
i?03 = R30 = R33 = 2uj 2 /c 2 . They coincide with the respective coefficients 
of the stress tensor: Too = T03 = T30 = T33 = 2uj 2 /c 2 . Therefore, once 
again, the Einstein equation is verified. The wave is free and we have R = 
and G = 0. 

Slightly more complicated is the case of a plane wave where divE is non- 
zero. This happens for instance when u = B2 = E\/c = f(x) sinco(t — z/c). 
As we know, the solution satisfies the equations (J3.1|) . ()3.2|) . (|3.3|) . but not 
the classical Maxwell equations. We suggest looking for a metric tensor of 
the form: 
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(13.16) 



where p is function of the variable £ = t — z/c. One has: \/—g = \fp\- The 
tensor is the same as in (j!3.7|) . Regarding the other tensors, we get: 
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62 



We must point out an extraordinary fact: in the new geometry, the 4- 
divergence of the electric field turns out to be zero. As a matter of fact, by 
noting that u/f and p do not depend on x, one has: 



PE 



1 d{^g F ok ) 



-g 



-1 d(\f P \F 



■01 \ 



dx k 

—c 
fp 2 



d{u/f) 







(13.19) 



\fp\ dx fp 2 dx 

Thus, the mass tensor is still vanishing. The Christoffel symbols are a little 
different from the ones in (|13.5|) (in particular is not zero), but the 
coefficients of the Ricci tensor are exactly equal to those given in (|13,6[) . 
Then, the equation (|13.11|) must be modified as follows: 



p"p 



P 2 X 



(13.20) 



thereby admitting the same solution p obtained in the case of the plane wave 
at uniform density. The laws of motion are the same as in (|13,13j) . They 
tell us that u and p shift at the speed of light along the z-axis. They do not 
specify however the function /, which must be assigned through the initial 
conditions. 

We are ready to illustrate the case of a spherical wave. With the 
same notation of sections 2 and 4, we set the coordinates in order to have: 
(xo, xi, X2, x%) = (ct,—r, —<fi,— 6). Let us assume that E = (0, cu, 0) and 
B = (0,0, u), where u = ^f(4>) sinw(t — r/c). We recall that, to avoid sin- 
gularities at the poles, the function / is not allowed to be constant. So, this 
is similar to the case of a variable-density plane wave. We also have that 
(V°, V 1 , V 2 , V 3 ) = (c, — c, 0, 0). So, let us begin by giving the metric tensor: 
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where p is a function of the variable £ = t — r/c. Note that here the case 
g22 = — 1 corresponds to the standard spherical system of coordinates. For 
the electromagnetic tensors we get: 



F ik 
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In order to evaluate i* 1 ^, we started from (|6.3|) . recalling that, by 13f) . 
one has (3>,A) = (— F(<p) sina;(t — r/c), — F{4>) sinu>(i — r/c), 0, 0), where 
F is a primitive of /. The metric in (|13.21|) is the same as the one we 
would have obtained if we had worked with a plane electromagnetic wave. 
The fact that we are in spherical coordinates is actually contained in the 
electromagnetic tensors (in which we find ru in place of u). Note that 
(Vb, V\, V2, V3) = (c, c, 0, 0), from which one obtains the relation F = 0, 
confirming that the wave is free. As far as energy is concerned, we get: 



/ {ru/pf) 2 (ru/pf) 2 
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Therefore, from the Einstein equation, we come to: 
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V X 
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(13.22) 



where we observe that the right-hand side only depends on the variable 
£ = t — r/c. The equation (|13.22f) once again gives the solution p = 
(vVx/u) sinu;(i- z/c). 

Finally, by differentiating T lk (for i = and i = 1), we find the Euler 
equation in spherical coordinates: 



1 du 
~c~di 



+ 



1 d{ru) 
r dr 







(13.23) 



The function u is the solution to (|13.23|) . after assuming the appropriate 
initial conditions. 

The results of this section, although only restricted to the analysis of 
free waves, bring to attention some important issues. Up to now, we have 
claimed that a good theory of electromagnetism was meaningful only by 
allowing divE to be different from zero. Here instead we find that pe = 0. 
In our opinion, what is happening can be explained as follows. The space- 
time "reacts" to the passage of a wave, by varying itself in syncronism, in 
order to make the 4-divergence of the electric field vanish. The perturbation 
of the geometry is however weak enough to maintain field curvature G equal 
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to zero. The classical divergence divE may instead attain arbitrary values. 
We ask ourselves if it is possible to set up an experiment showing that, 
at some point and at a certain time, one has divE 7^ 0. Perhaps, this is 
not possible since, due to the modification of the metric, the instruments 
are unavoidably affected by the deformation of time and distances (with 
respect to the Euclidean reference frame). Therefore, in place of divE, we 
could end up measuring pE- But the last quantity is always zero (a least for 
free waves). As a consequence, we conclude that some divergence vanishes, 
although is not the classical, but the relativistic one. According to this new 
interpretation of the facts, in some sense the Maxwell theory was correct. 

There is another point that needs to be clarified. The problem is why 
the geometry changes depending on the electric field, and not the magnetic 
field, expecially after we said that for free waves the two fields have the 
same role. Firstly, we note that, in all examples studied in this section, 
the condition divB = was always fulfilled. In addition, if similarly to pe> 
we define p~B, we discover that this new quantity is also zero (see section 
14). If we imagine the wave like a fluid in motion, then this condition says 
that there is no flow of some "magnetic density of matter". In truth, it 
is reasonable to assume that a sole electromagnetic fluid exists (not two, a 
separate electrical one and a magnetic one). As will become clear in the 
next section, where we analyze the case divB / 0, such a fluid turns out 
to pulsate along a specific tangential direction (in principle, not necessarily 
corresponding to that of the electric field). Exchanging cause with effect, in 
section 14 we will support the following statement: from the behavior of the 
natural events, we are inclined to name the direction of the electric field as 
being that identified by the transversal oscillations of the fluid in motion. 

Far more complicated phenomena show up, when we suppose that the 
waves are no longer free (thus, G/0). In this context, the real gravitational 
fields come into life. We do not have any specific examples to discuss, due to 
the difficulty of the problems involved. Some hints will be given in section 
15. 



14 The divergence of the magnetic field 

In the previous sections, some situations were discussed under the hypoth- 
esis divB = 0. Although our equations now have a general validity, the 
assumption is necessary, for instance when introducing the potentials A 
and <E>. Regarding this condition, we would like to add further comments 
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in this section. It is standard to introduce a transformation that exchanges 
the role of the electric and magnetic fields. This can be done through the 
pseudo-tensor: 

{0 when at least two indices are equal 
1 if the indices form an even permutation (14.1) 
— 1 if the indices form an odd permutation 

The parity of the permutations is counted starting from the set: {0, 1, 2, 3}. 
Then, we define: 

£ ''" = ^m^j^i^k^mjik = ^rajik (1^.2) 

We may now introduce the duals of the tensors (|6.4j) and (|6.5|) in the fol- 
lowing way: 

F-mj — "2^-mjikP F ^ — d m Cj F m j — ^ ^ik (l^ - ^) 

Therefore, we obtain for example: 

77 7723 77 773I 77 J?12 77 7701 771 Z?02 77 7703 

i*01 — f 1*02 — r -^03 — # ^23 — r #31 — # #12 — # 

The original tensors and their duals have the same structure, with the dif- 
ference that E replaces — cB and cB replaces E. 

In a similar way, the dual of the anti-symmetric rank-three tensor Fjik 
(defined in (|11.13[) ) is given by: 

pm = _X e mjik F .. k (144) 

Hence, up to even permutations of the lower indices, one has: 

F~0 77i 77 1 77 772 77 773 77 

— ^123 # — —-^023 # — #013 # — — #0\2 

In general coordinates, it is customary to define: 

r-mjik / ~~ „mrn'jj'ii'kk' 1 

G — y-9 9 9 9 9 tm'j'i'k' — - ^j—^ tmjik 

^mjik = \J 9 ^mjik (14.5) 
So that the duals in (|14,3|) and in (|14.4|) are generalized as follows: 

p ■ — I (= ... F ik F m i — —I c m i ik F m — —1 ^mjik F fi \ 

1 raj — 2 ^nijik 1 1 — g — 6 1 Jik ) 

where F^ is provided in (|6.4|) and F* can be found in H12.11|) . 
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Then, the following relation is known (see [7], p. 134): 

1 d{^~g F m i) ~ m 



d Xj 



;i4.7) 



where we supposed that F m is the dual of the cyclic derivative Fj^ of the 
tensor F^ (of which F m:> is the dual). 

Passing to the duals, the equation (jll.27f) becomes: V°F m = V m F°. 
Therefore, by (|14,7|) we get: 
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m 



0,1,2,3 (1U 



which is the exact counterpart of l|11.24|) . The equation 1)14. 8() represents, in 
a general coordinates system, the equation (|9.5|) . that is equivalent to (|9,4|) . 
after taking E in place of — cB and cB in place of E. From (|14.8|) . we can 
recover the continuity equation: 



1 diV^gfmV* 



-9 



with pb 



1 d{^9 



-a 



dx k 



(14.9) 



It is worth noting that p& has the same dimensions of /9e- For example, 
according to (|14.6|) . the dual of Fn- in (|13.7|) is: 
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(14.10) 



Since we supposed that u does not depend on y, we obtain p-& = 0. Based 
on the metric given by (|13.16|) (where \J —g = \fp\), we just checked that, 
together with the condition divB = 0, the 4-divergence of the magnetic field 
also vanishes. 

Going back to the equation ()12.24|) . this time we cannot assume that 
Fmjk = 0. On the other hand, we can use (fT4~71) and (fH~%]) with V° = c, to 

get: 



-n im F uF jk 

— 9 rimV 



g im F l F lm 



1 3(^=5 FV) 



-r gimV 



dxj 



— F im Vrr, 



Fir. 
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The first passage follows on from a direct counting of the permutations of 
the indices, thanks to the definitions provided in Q14.6JI . Substituting in 
(|12.24[) . one finally gets: 

V k W k = - c ( PE F im V m + pvF im V m ) (14.12) 

Let us observe that, for i = 0, we have F 0m V m = (due to (pTTnjl V while 
pnF 0m V m recalls the product -c 2 (B- V)divB. Thus, the first line of (114.12)) 
turns out to be equivalent to 1)10.4)1 . 

The equation (|14.12|) is the generalization of (|12.6|) with N = and 
M = 0. In spite of its elegance, it is not very convincing, since it involves 
two mass densities (see also the comments at the end of section 13). Let us 
try to explain what is happening. Without going into technical detail, we 

may make some remarks. We first note that F{ k = —Fi k , which means that, 
after applying the dual twice, one gets the opposite of the original tensor. 
Then, for any real A, we consider the two tensors: 



V* 2 + (1 - A) 2 



XF lk + (1 - X)F 



ik 



fik 



XF l « - (1 - \)F ik (14.13) 

v/A 2 + (1 - A) 2 L 1 ; lk \ K ' 

where the second one is the dual of the first one. As in 1)12.11)) we have: 
jrik _ g im g kl F ml . Moreover, we can check that the tensor t7 ifc in ()12.1())) 
does not change if in place of F^ and F lk we take J~i k and J- ik , respectively. 
Therefore, the electromagnetic stress tensor does not depend on A, even if 
this parameter varies in space and time. Actually, we already observed in 
section 12 that the energy tensor does not recognize the polarization of the 
electromagnetic field. 

At this point, we can introduce the two new densities (see also ()12.7)) ): 
pe = V k T ok and p B = V k T ok , where £ = (A E + (1 - X)cB)/ v / X 2 + (1 - X)' 2 
and B = (AcB - (1 - A)E)/^A 2 + (1 - A) 2 . 

So, another equivalent way to write equation 1)14.12)1 is: 

V k U lk = - c ( P£ F m V m + p B P m V m ) (14.14) 

For A = 1 the two versions are actually the same. Now, by letting A to 
vary, suppose that it is possible to modify the polarization of the fields £ 
and B at each point, in order to get p B = 0. In this way, we are left with 
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a single density pg, which is the one to be used in constructing the mass 
tensor M ik = p £ ViV k . 

Let us restate the situation in brief. Every non-trivial electromagnetic 
wave presents regions where the classical divergence of any of the two fields 
is non-zero. The eletromagnetic energy tensor does not distinguish between 
the two types of fields (electric or magnetic). In the end, what matters is the 
intensity of the wave and the modality of propagation of its fronts, without 
paying attention to the way each front has been parametrized. We can 
associate a fluid in motion at the speed of light with the wave. Independently 
of the actual orientation of the fields E and B, we can locally build two other 
fields £ and B, so that the first one oscillates together with the fluid and 
the second one satisfies divB = 0. This fictitious change of polarization 
has no influence on the electromagnetic energy tensor. The 4-divergence of 
£, when different from zero, represents the mass density of the fluid and 
it is used to construct the mass tensor. This last tensor is added to the 
electromagnetic energy one, to form the global energy tensor which is on 
the right-hand side of the Einstein equation. In principle, the fields £ and B 
are not directly associated with E and B. However, in the natural evolution 
of electromagnetic phenomena, the two entities usually coincide. 

All the examples analyzed in the previous section satisfy pjs = and 
pB = 0, hence, they were already well suited to the case A = 1, correspond- 
ing to £ = E and B = cB. In particular, the case in spherical coordinates 
simulates the real behavior of a wave generated by an infinitesimal electric 
dipole ascillating in a vertical direction. Somehow, the dipole imparts me- 
chanical oscillations to the fluid, in the same direction as the electric field. 
Formally, we can now exchange the role of the fields E and cB, by polarizing 
the spherical wave by 90 degrees. In this new situation, we have divE = 0, 
pE = and divB ^ 0. By choosing A = 0, we realize the condition ps = 
and the fictitious fields £ and B turn out to be anti-ruotated by 90 degrees. 
Therefore, there is no longer coincidence of £ and B with the corresponding 
E and cB. Nevertheless, a spherical wave having the second kind of polar- 
ization is difficult to observe in nature, since it should correspond to the one 
generated by an infinitesimal magnetic monopole. 

It is certainly true that our equations are not capable of recognising the 
polarization of free waves. This is a property that comes with the initial 
conditions. However, free waves are created by some causes inherent to nat- 
ural events, which have a strong influence in determining polarization. The 
problem resides at the origin, for example in the non existence of magnetic 
monopoles (we will have a short discussion about this in section 15). Recall 
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that, in equation Q9.6JI . the electric and magnetic fields cannot be inter- 
changed. Certainly, this equation influences the creation of a spherical wave 
through the mechanical oscillations of an electric charge. The conclusion 
is that, at least for free waves, we can expect A = 1, which implies that 
the direction of transversal propagation of the fluid is in accordance with 
that of the electric field. More precisely, this can be taken as a definition of 
electric field. Suppose that an external mechanical perturbation is applied 
to a free wave having £ = E, in a direction not alligned with that of field 
E, in such a way the direction of £ changes. Then we may think that the 
wave reacts by varying its polarization (see section 7, 8 and 9) in order to 
correct its posture, bringing field E to once again coincide with field 8. In 
other words, the electric field turns out to be identified with the one that 
follows the transversal oscillations of the fluid, and such a definition matches 
reality. 



15 Other developments and conclusions 

We start by recalling the primary results obtained by the paper. In section 
9, we introduced the following equations: 

i9E 

— = c 2 curlB - (divE)V (15.1) 

dt v ' y ' 

— = -curlE - (divB)V (15.2) 
DV 

— = ^(E + V x B) (15.3) 

where E is the electric field, B is the magnetic field and V is a velocity field 
satisfying: 

|V| = c (15.4) 

The constant fi is a charge divided by a mass, and c is the speed of light. 

Then, in sections 11 and 12, we wrote the equations in covariant form. 
In the same order they appear above, we have, for i = 0, 1, 2, 3: 

(V k F ik )V° = (V fe F ofe )F* (15.5) 

(V m F m )V° = (V m F 0m )V l (15.6) 

—- + Ti k V*V k = -!=■ F m V m (15.7) 
Vt J c 
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where V° = c, Fi k is the electromagnetic tensor and F^ m its dual. In the 
general system of coordinates, the normalizing condition takes the form: 



The metric tensor is not given, but has to be determined through the 
Einstein equation: 



= ^(-M mj F im F kj + l mk F mj F m i + ViV k V m F 0m ) (15.9) 



where on the right-hand side we find a suitable energy tensor, obtainable 
with the rules provided in section 12 (in the construction of pe = V m -F 0m 
remember to take into account the warnings at the end of section 14) . Such a 
coupling corresponds to a quite complex system, able to describe space-time 
geometry in conjunction with electromagnetic phenomena. 

Our set of equations contains the embryo of some of the main laws of 
Physics. In (|15.1|) and (|15.2|) . we recognise the equations of electromag- 
netism, more or less with the same structure as the Maxwell equations. We 
discovered that, when V is irrotational, then (|15.4|) is the eikonal equation, 
so that the Huygens principle is also latent. On the right-hand side of (j!5.1|) 
we partly recognize the Ampere law. The equation 1)15. 3|) expresses the 
Lorentz law, anticipating the Newton law in the form of momentum equa- 
tion for the dynamics of fluids. In fact, we claimed that the light rays can 
be assimilated to stream-lines of a certain fluid of density pe- Moreover, we 
know that a continuity equation holds for pe- 

Throughout the paper, we assumed we were in a universe that we could 
call "pre-Coulombian" . As a matter of fact, we developed a theory of elec- 
tromagnetism without introducing any charges, and we spoke about fluids 
without having any masses. The only elements at our disposal were the 
fields. Here comes the big question: can we now build matter from these 
fields respecting the rules that we wrote? In other terms: can an elementary 
particle be "solution" to our set of equations? 

A particle is a quite complicated thing. It has charge, magnetic mo- 
mentum, spin, mass. It evolves and interacts with other particles according 
to the rules of quantum mechanics. Can we contain all these factors in 
a solution localized in space? This problem was mentioned in section 10 
where we discussed possible solutions, consisting of a stable system of two 
rotating solitons. Although the framework is still incomplete we collected 
some pieces of evidence, whose details will shortly be discussed below, that 



g im V l V m = 
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support the possibility of creating particles from fields. We recall that other 
authors, through a qualitative analysis, followed a similar idea of building 
electrons from photons (see for example |13j and the references therein). 

We can give a rough idea of how a "particle solution" looks by examining 
figure 7, that shows, projected on a plane, the rotation of the fronts around 
an axis. From a qualitative viewpoint, field E oscillates radially, but, in the 
average, mainly pointing inward (or outward). This creates the polarity of 
the electrical charge. Field B is orthogonal to the page. The rays form closed 
orbits and their vector curvature G points toward the center, producing a 
non-vanishing gravitational field. If the sign of E is changed, then G again 
points toward the inside (gravity has only one polarity). The displacement 
of field V matches the idea that something is "spinning" , and the associated 
electromagnetic fluid corresponds to a kind of vortex. 

Still referring to figure 7, let us suppose that the particle is an electron. 
Then E should be directed toward the center and, using the standard vector 
product x, B points downwards. Nevertheless, a negative charge rotating 
clockwise produces a spin angular momentum pointing downward and a 
magnetic field pointing upward, which is in contrast to what previously 
found. As we remarked in section 8, this happens because we do not use 
the suitable vector product x . In fact, the correct one is the left-handed 
one. Since the magnetic dipole moment is independent of the sign of x , the 
change of parity now confirms that B points upward. If we want to maintain 
the same set of equations, we can solve the problem just by changing the 
sign of the electric field, so that the electron has a chance of existing only 
if the electric field vectors point outward. We can still call this particle an 
"electron" and give a negative sign to it, but we have to comply with the 
new rule stating that currents flow from a negative pole to a positive one. 

Nevertheless, one can see that such a situation is still not compatible 
with equations ((To^ - (fTo^|) - lfrO)) - ((TO|) . One of the reasons is that the 
outer orbits of the light rays are longer than the inner orbits, and this does 
not match the condition (|15.4|) . telling us that the information propagates 
at constant speed. In order to have chances of finding solutions to the form 
described above, the use of the general relativity framework is unavoidable. 
The modification of space-time geometry allows for the preservation of the 
momentum of inertia (a typical mechanical concept), providing the "glue" 
that keeps the particle together. The rotating wave follows the geodesies in 
the new metric. At the same time, the curvature of such geodesies has to be 
compatible, through (|15.7|) . with the electromagnetic setting. The geometry 
alters the relation between space and time in such a way that the rays, always 
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travelling at speed c, can accomplish paths of different length in the same 
amount of time. This recalls the problem of the rigid rotating disk in general 
relativity. It is clear that the particle solution (if it exists) involves the use 
of the whole set of equations. Therefore, its determination, even from the 
point of view of numerical computations, is a demanding problem. Finally, 
by heuristic arguments, one can recognise that a similar solution, where the 
magnetic field is exchanged with the electric one, should be forbidden by 
equation (|15.7|) . This would imply the impossibility of building magnetic 
monopoles. 

We finish the paper with some further speculations, not having enough 
theoretical background. One positive aspect is that particle solutions are 
expected to be extremely stable (an electron is quite a difficult object to 
destroy). Another aspect is that they are in some sense "unique" (there is 
only one type of electron or proton), and this property raises other ques- 
tions. The equations and (|9.2fl are "scalable", by meaning that we 
can multiply the fields of a free wave by a constant, oncemore obtaining a 
solution. Thus, free waves may be of any size and intensity. But, if we take 
into account constrained waves, then this property is no longer true, since 
(|15.4|) is not a scalable equation. Together with c, \x and x m <113.1|) . another 
constant is hidden in the set of equations (|15,5j) - (fT5.6j) - (|15.T|) - ()15.8j) . which 
is related to some "magnitude" of the geometry. This result does not penal- 
ize our theory. Actually it may give more strength to it. As a matter of fact, 
we cannot have electrons of any size! We have no elements for quantifying 
the values of the various parameters, unless we find the particle solution 
explicitly. 

The last issue we discuss is the convenience of setting up experiments 
validating our theory. The problem is left to the experts. However, we 
think that many convincing arguments, also based on a multitude of prac- 
tical observations, have already been collected, showing that our model is 
adequate. The real breakthrough would be in predicting the realization of 
an electromagnetic device, capable of producing gravitational field. 
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